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PREFACE 

•The  work  in  this  report  was  performed  under  Contract 
No.  F40600-74-C-0013  for  the  Arnold  Engineering  Develop¬ 
ment  center,  Arnold  Air  Force  Station,  Tennessee,  37389, 
by  Block  Engineering,  Inc.,  19  Blackstone  Street,  Cam¬ 
bridge,  Massachusetts,  02139.  This  report  includes  all 
work  performed  during  the  contract  period  22  April  1974  to 
April  1975  under  Program  Element  65807F.  The  report  was 
submitted  on  1  October  1975.  The  work  was  performed  under 
the  direction  of  Air  Force  Technical  Manager,  Kenneth  B. 
Harwood,  Captain  (CF) ,  Directorate  of  Technology,  Arnold 
Air  Force  Station,  Tennessee. 

The  authors  acknowledge  the  generous  help  of  Dr. 
Stewart  Hemple  of  Block  Engineering,  Inc.  and  Professor 
Ellis  R.  Lippincott  (deceased)  of  the  University  of  Mary¬ 
land  in  the  theoretical  portions  of  this  report. 

The  reproducibles  used  in  the  publication  of  this 
report  were  supplied  by  the  authors. 
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SECTION  I 
INTRODUCTION 


1.0  GENERAL 

This  is  the  Final  Report  for  Contract  No.  F40600-74-C-0013, 
with  Arnold  Engineering  Development  Center,  Arnold  Air  Force 
station,  Tennessee.  This  Contract  contains  both  experimental 
measurements  of  the  Raman  Stokes  scattering  cross-sections  of 
water  vapor  as  a  function  of  temperature,  and  a  theoretical 
calculation  of  the  spectral  profiles  for  the  asymmetric  top 
model  of  the  water  vapor  molecule  as  a  function  of  temperature. 

Laser  Raman  scattering  diagnostic  techniques  are  being 
developed  at  AEDC  for  the  measurement  of  concentration,  and 
vibrational,  rotational  temperatures.  This  is  to  be  applied 
to  single  and  multiple  gaseous  species  in  both  single  and 
multiphase  flow  fields  for  flow  diagnostic  test  applications. 
These  fields  occur  in  the  exhausts  of  jet  engines  and  rockets 
and  in  combustion  chambers.  Water  vapor  is  present  to  some 
degree  as  a  combustion  product  in  all  these  applications. 

One  phase  of  this  program  entailed  the  development  of  two 
burners  which  produced  H2/O2  and  H2/Air  flames.  These  burners 
produced  steady  temperature  and  flow  rate  flames.  Thus,  water 
molecules  in  thermal  equilibrium  were  obtained  over  a  suffi¬ 
ciently  large  area  to  allow  their  examination  by  laser  Raman 
spectroscopy. 

1.1  ORGANIZATION  OF  THE  REPORT 

This  report  is  divided  into  two  major  parts.  The  first. 
Sections  II  and  III,  describes  the  experimental  apparatus, 
procedure  and  results.  Section  II  contains  information  on  the 
spectrometer,  burner  design,  flow  meters,  thermocouple  data 
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and  the  experimental  procedure.  Section  III  presents  the 
spectral  data  and  a  discussion  of  the  experimental  error. 

The  second  part.  Sections  IV  and  V,  discusses  the  theoreti¬ 
cal  considerations  and  calculations.  The  results  of  these 
calculations  are  tabulated  in  Section  V.  The  computer  program 
used  in  the  theoretical  calculations  can  be  found  in  the 
Appendix. 
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SECTION  II 

•  INSTRUMENTATION  AND  EXPERIMENTAL  METHOD 

2.0  GENERAL 

This  section  describes  the  Raman  spectrometer  used  in 
acquiring  all  the  spectral  measurements  presented  in  this 
report.  It  describes  the  design  and  development  of  the  burners 
and  the  calibration  of  the  flow  meters.  Discussions  are  devoted 
to  temperature  measurements  by  thermocouple  for  determination 
of  temperature  gradients  in  the  flame  and  use  of  the  Raman 
profile  of  the  N2  spectrum  for  obtaining  flame  temperatures. 

The  experimental  procedure  is  described  in  detail  for  the 
manner  in  which  all  measurements  were  made. 

2 . 1  RAMAN  SPECTROMETER 

The  initial  spectral  observations  are  confined  to  the 

O 

Stokes  bands  arising  from  4880  A  incident  radiation  from  an 
argon  ion  laser  (Spectra  Physics  Model  165)  operated  for  all 
data  at  1.2  watts.  The  scattered  light  was  analyzed  by  a 

O 

double  monochromator  (Jarrell  Ash  Model  25-100)  with  5000  A 
blazed  gratings.  The  detector  was  a  cooled  ITT-FW130  photo¬ 
multiplier  operated  in  the  pulse- counting  mode  with  dark 
count  levels  of . about  5  counts  per  sec  for  this  work. 

The  overall  experimental  arrangement  was  designed  to 
have  the  laser  beam  travelling  parallel  to  the  entrance  slit 
(that  is,  vertically)  such  that  the  Raman  scattered  light  is 
collected  90®  to  the  incident  light  direction.  The  Raman  scat¬ 
tered  light  was  collected  by  a  lens  with  focal  length  100  mm 
(F  number,  8)  .  The  width  of  the  laser  beam  was  eOaout  100  ym. 

One  to  one  magnification  of  an  image  of  1.5  cm  in  height  from 
the  scattering  zone  area  was  accepted  in  the  entrance  slit 
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of  the  monochromator.  The  slits  were  set  at  350  ym  which 
produced  a  spectral  slit  width  of  about  8.0  cm  ^  for  the  wave- 
Tength  of  the  scattered  light  considered. 

In  following  our  regular  procedure,  the  Jarrell- Ash  Model 
25-100  Raman  monochromator  was  aligned  for  optimum  throughput. 
Initial  calibration  of  the  spectrometer  was  carried  out  by  a  • 
representative  of  Jarrell-Ash.  Prior  to  starting  the  collection 
of  experimental  data  we  rechecked  the  calibration  of  the 
monochromator.  The  wavelength  of  light  used  and  the  sources 
from  which  they  were  obtained  are  listed  in  Table  1. 

TABLE  1 

MONOCHROMATOR  WAVELENGTH  CALIBRATION 


Source 

O 

X,  A 

Grating 

Order 

cm 

Hg  lamp 

4358.3 

I 

22944.5 

5460.7 

I 

18312.5 

4046.6 

II 

12356.2 

Argon  laser 

4579. 

I 

,21839. 

4880. 

I 

20492. 

5145. 

I 

19436. 

Figure  1  is  a 

plot  obtained 

from  this  data. 

Wavenumber 

difference  (difference  between  true  wavenumber  and  absolute 
wavenumber  reading  on  monochromator)  is  plotted  against  true 
wavenumber  of  the  source. 

2.2  BURNER  DESCRIPTION 

2.2.1  Porous  Plugburner 

The  flames  studied  in  the  1000  to  1500  K  range  were 
produced  by  a  water-cooled  sintered  brass  porous  50  ym) 
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plug  burner  (diameter,  2.5  cm).  It  was  operated  horizontally 
and  burned  into  another  water-cooled  porous  plug  placed  about 
2  cm  away  which,  in  turn,  was  connected  to  a  rough  vacuum  line. 
In  this  fashion,  a  stable  horizontal  flame  at  atmospheric  pres¬ 
sure  was  produced.  There  are  several  advantages  to  this  type 
of  burner;  (1)  the  laminar  flame  is  fairly  uniform  over  a 
rather  large  area  offering  an  ideal  scattering  test  zone,  (2) 
the  outer  layers  of  the  flame  have  been  shown  to  serve  as  a 
shield  to  inhibit  atmospheric  contamination  of  the  inner  flame 
which  will  be  studied.  The  burner  is  shown  in  Figure  2.  Con¬ 
tained  in  Figure  3  is  a  top  view  of  the  burner  setup. 

An  excellent  description  and  study  of  this  burner  design 
can  be  found  in  the  literature  on  combustion.  The  main 

advantage  of  this  burner  is  that  it  is  possible  to  cool  a  flame 
to  a  lower  temperature  than  that  which  would  result  from  adia¬ 
batic  combustion.  The  resulting  flames  are  very  stable  and 
extremely  flat  -  independent  of  the  size  of  the  burner.  The 
previous  studies of  this  burner  have  also  demonstrated  that 
flame  temperature  can  be  varied  approximately  200  K  by  merely 
changing  the  pre-flame  gas  flow  rate  and  keeping  the  gas  con¬ 
stituents  constant.  See  Table  2  for  some  results  on  H2/  O^t 
and  flames. 
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TABLE  2 

EXPERIMENTAL  DATA  SHOWING  MEASURED  FLAME  VELOCITY  VERSUS 
FLAI/IE  TEMPERATURE  DATA  FOR  CONSTANT  AND  PRE-FLAME  COMPOSITION* ** 

Flame  1,  2H2 ;  5N2:102 

Velocity*  (mA/s) _ Temperature  (K) 


59.0 

1550 

98.5 

1667 

173.0 

1818 

..2H2;  5N2;  IO2 

Velocity  (ml/s) 

Temperature  (K) 

49.2 

1409 

108.0 

1540 

197.0 

1667 

*  Velocities  are  evaluated  for  the  total  pre-flame  compositions 
at  the  stated  molar  ratio  at  25°C  and  1  atmos.  pressure. 

**  Data  collected  by  W.  E.  Kaskan  (Ref.  1) . 


17 


AEDC-TR-75-1  50 


2.2.2  Needle  Burner 

•Due  to  the  characteristics  of  the  flame  produced  by  the 
porous  plug  burner,  it  was  necessary  to  construct  another 
bxirner  in  order  to  obtain  reasonably  good  data  at  temperatures 

I 

higher  than  1500  K.  The  water-cooled,  metal,  porous  plate 

burner  produces  a  higher  temperature  only  in  the  flames  tip 

and  only  over  a  very  small  region  which  is  not  satisfactory 

for  obtaining  good  Raman  spectra.  The  alternate  burner  is 

constructed  from  400  stainless-steel  hypodermic  tubes  tightly 

soldered  together.  The  internal  diameter  of  the  tubes  is 

0.045  cm.  This  gives  primary  cones  not  more  than  1/2  mm  in 

height,  and  an  excellent  laminar  flame.  The  experimental 

flame  is  about  2  cm  in  diameter.  The  specific  design  is  after 

(2) 

that  of  Padley  and  Sugden  .  The  burner  is  water  cooled  and 
premixes  the  gases.  The  burner  easily  obtains  temperatures 
between  1500  K  and  2500  K.  Experimental  studies  have  shown 
that  the  flame's  central  portion  is  well  shielded  up  to  10  cm 
away  from  the  burner  head.  This  burner  is  shown  schematically 
in  Figure  4 . 

2.3  DESCRIPTION  OF  THERMOCOUPLE 

% 

The  thermocouple  was  prepared  from  platinum  and  platinum-13% 
rhodium  wire  of  0.0005  inches  (0.0127  mm)  in  diameter.  The 
thermocouple  junction  was  butt  welded.  The  wire  is  contained  in 
ceramic  insulating  tubes.  The  junction  is  supported  above  the 
insulating  tube  such  that  only  the  thermocouple  wire  is  in  the 
flame  area.  Low  velocity  and  low  temperature  H2-air  flames  are 
known  to  produce  very  little,  if  any,  catalytic  heating  at  the 
platinum  thermocouple's  surface ,  so  it  will  not  be  necessary 
to  coat  the  thermocouple  with  fixed  quartz  or  ceramic  materials 
to  eliminate  this  effect. 
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TUBE  BURNER 


STAINLESS  STEEL  TUBES 


Cu 


FIGURE  4 
TUBE  BURNER 
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Kaskan^^^  has  shown  that  corrections  for  radiation 
cooling  caused  by  the  thermocouple  wire  are  significant  and 
may  be  calculated  from 


AT 


1.25eaT 

w 


rad 


0.25 

_n 

pv 


(2-1) 


where 


e  is  the  thermocouple  emissivity, 
a  is  the  Stefcin-Boltzmann  constant, 
is  the  wire  temperature, 

D  is  the  wire  diameter, 

X  is  the  thermal  conductivity  of  the  gas  at  T^, 
n  is  the  viscosity  of  gas  at  T  ,  and 

Vv 

2 

pv  is  the  mass  fl\ix,  gm/sec/cm  . 


Estimates  of  AT  ,  were  obtained  using  the  values  of  the  para- 

(1) 

meters  in  equation  (2-1)  as  given  by  Kaskan'  '  for  temperatures 

between  1000  K  and  1800  K.  The  results,  however  crude,  are 

(3) 

comparable  with  those  of  Lapp'  and  are  given  below. 


Measured  Temperature 


Correction  Due  to  AT 


rad 


1000 

K 

18 

1200 

K 

30 

1400 

K 

40 

1600 

K 

50 

1800 

K 

65 
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2.4  DETERMINATION  OF  TEMPERATURE  FROM  THE  CALCULATED  RAMAN 

SPECTRA  OF  N2 

In  Figure  5  is  a  series  of  calculated  Raman  spectra 
of  the  vibrational  Stokes  scattering  of  N2  at  various  tempera¬ 
tures  showing  the  intensity  dependence  on  temperature  of  the 
first  vibrational  hot  band.  The  calculations  are  similar  to 

(3) 

those  previously  described  in  the  literature  .  Plotted  in 
Figure  6  is  the  ratio  of  the  areas  of  the  N2  fundamental  and 
its  hot  band  versus  temperature.  By  using  the  observed  relative 
intensities  of  N2/  the  vibrational  temperatures  of  N2  can  be 
found  from  the  curve.  As  will  be  observed,  these  temperatures 
agree  remarkably  well  with  those  of  the  corrected  thermocouple 
measurements  and  they  probably  have  the  same  accuracy  of  ±30 
to  40  K.  Use  of  a  thermocouple  is  definitely  required  in  order 
to  insure  that  the  Raman  scattering  zones  utilized  contain  only 
small  temperature  gradients. 

2.5  DESCRIPTION  OF  FLOWMETERS 

The  flowmeters  used  in  this  study  are  commercially  avail¬ 
able  Roger  Gilmont  Flowmeters.  The  models  used  consisted  of 
two  size  number  3  and  one  size  number  2.  Size  number  3  has  a 
useful  range  of  200  to  1,200  ml/min  of  air.  The  number  3  flow¬ 
meters  were  used  to  monitor  the.  H2  and  N2  flowrates.  The  number 
2  flowmeter  was  used  for  the  O2  flowrate. 

In  Figure  7  is  shown  a  comparison  of  a  volume  displace¬ 
ment  calibration  curve  and  a  calculated  curve.  The  calculations 
are  based  on  a  method  supplied  by  the  manufacturer.  The  experi¬ 
mental  and  calculated  curves  agree  very  well. 

In  Figure  8  are  volume  displacement  calibration  curves  com¬ 
paring  two  different  flowmeters  of  the  same  type.  Agreement 
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FIGURE  5 

CALCULATED  SPECTRUM  AT  VARIOUS  TEMPERATURES  (DEGREES  K) 
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FIGURE  7 

CAL IBRATION  OF  FLOWMETERS 
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here  is  very  good,  indicating  reproducibility  in  the  manufacturing 
of  the  flowmeters.  Water  was  the  liquid  used  in  the  volume 
displacement  experiments  and  the  appropriate  corrections  for 
the  vapor  pressure  of  H2O  and  the  solubility  of  the  gases  in 
water  have  been  applied  to  the  results. 

Because  of  the  good  agreement  between  the  theory  and 
experimental  results  and  in  order  to  avoid  experimental  errors , 
calculated  calibration  curves  were  used  for  determining  flow 
rates.  Presented  in  Figures  9,  10  and  11  are  experimental 
calibration  curves  for  H2,  N2  and  O2  respectively. 

2.6  EXPERIMENTAL  PROCEDURE 

The  first  step  in  obtaining  any  spectral  data  consisted 
of  accurately  aligning  the  laser  beam  and  burner  with  respect 
to  the  slits.  Also  included  in  this  process  is  the  accurate 
mapping  of  the  flame  as  a  function  of  temperature.  Since  the 
thermocouple  was  used  an  an  indicator  of  the  foregoing  pre¬ 
requisites,  it  was  necessary  to  mount  it  on  a  precisely  cali¬ 
brated  x-y-z  translator.  This,  in  turn,  was  mounted  on  the 
optical  rail.  A  low  power  laser  beam  was  then  turned  on  and 
the  usual  alignment  procedure  was  carried  out  to  position 
the  beam  parallel  to  the  slits.  The  center  of  the  thermocouple 
wire  was  then  aligned  with  the  laser  beam,  so  that  a  vertical 
path  of  the  thermocouple  was  parallel  to  the  entrance  slit. 

The  laser  beam  was  then  blocked  off  and  the  burner 
ignited.  The  flame  was  adjusted  to  the  approximate  temperature 
to  be  investigated.  Vertical  measurements  of  temperature  were 
taken  about  every  2  mm  to  define  the  center  of  the  flame. 

Next,  by  manipulation  of  the  height  of  the  burner  and 
repetitive  temperature  measurements ,  the  center  of  the  flame 
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was  aligned  with  the  center  of  the  slit.  'Parallelism  and 
alignment  was  again  rechecked  by  use  of  a  low  power  laser 
beam. 

The  burner  was  then  turned  off,  the  thermocouple  lowered 
out  of  the  optical  path  and  the  spectra  of  the  nitrogen  and 
water  vapor  content  of  the  laboratory  atmosphere  were  obtained. 
This  was  carried  out  primarily  as  a  final  check  on  instriunent 
throughput.  Typical  spectra  of  this  type  are  shown  in  Figures 
12  and  13  respectively.  The  radiative  background  for  the 
spectra  is  about  20  cps. 

Without  disturbing  anything,  the  flame  was  then  reignited 
and  adjusted  to  the  appropriate  temperature.  The  parallelism 
and  alignment  of  the  thermocouple  with  respect  to  the  slits 
and  laser  beam  was  checked  again.  The  flame  was  then  mapped 
by  measuring  the  temperature  at  1  mm  intervals ,  for  a  distance 
of  5  mm,  on  each  side  of  the  center  of  the  flame.  The  slit 
heights  were  then  closed  to  a  height  of  10  mm.  A  low  power 
laser  beam  was  allowed  once  more  to  shine  on  the  thermocouple 
and  the  image  followed  on  the  slits  to  insure  that  the 
temperature  limits  to  be  investigated  fell  within  the  allowed 
slit  height,  and  that  those  not  desired  were  blocked  off. 

A  typical  example  for  further  discussion  would  be  that 
of  the  1500  K  flame.'  Here,  the  measured  temperature  of  the 
flame  in  the  area  where  Raman  scattered  light  was  collected 
varied  from  a  maximum  of  1550  K  at  the  center  to  as  low  as 
1425  K  at  the  outer  edges  of  the  measured  zone.  Several 
reasons  can  exist  for  these  variations.  Cold  air  is  possibly 
mixing  with  the  less  dense  hot  flame  gases  to  cause  the  flame 
to  break  down  and/or  there  is  a  large  inhomogenous  flame 
produced  from  irregular  surface  pores  causing  an  uneven  flow 
of  the  gases.  Before  obtaining  accurate  Raman  cross-sections 
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FIGURE  12 

ROOM  TEMPERATURE  RAMAN  SPECTRUM  OF 
IN  LABORATORY  ATMOSPHERE 
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of  water  vapor  in  various  flames ,  it  is  necessary  to  completely 

characterize  the  flames  across  the  burner  surface.  Previous 
12) 

studies  '  of  flames  produced  with  similar  burners  have  shown 
that  it  was  necessary,  in  order  to  obtain  a  homogenous  tempera¬ 
ture  region  in  the  flame,  to  collect  scattered  light  from  only 
a  5  mm  section  of  the  flame.  Also,  increasing  the  mass  flow 
of  the  flame  will  produce  a  more  stable  flame  with  less  atmos¬ 
pheric  influence.  Previous  studies have  also  shown  that 
by  placing  a  steel  screen  immediately  cdsove  the  burner  a  much 
flatter  and  more  homogenous  flame  can  be  produced. 

It  is  quite  evident  in  the  spectral  data  obtained,  even 

at  room  temperature,  that  asymmetry  exists  in  the  vibrational 

bands  of  water  vapor.  This  asymmetry  is  observed  to  increase 

sharply  in  the  1500  K  flame.  The  half -band  width  of  the  water 

band  at  room  temperature  approximately  doubles  at  the  higher 

temperature  and  the  maximum  peak  position  decreases  in  energy 

from  'V  3652  cm”^  to  3646  cm  A  shoulder  on  the  low  energy 

side  of  the  band,  (an  15  cm  ^  shift)  also  is  observed  to 

appear  at  the  high  temperature.  The  low  energy  asymmetry  is 

produced  by  the  fact  that  all  the  AJ  =  0  rotational  lines 

corresponding  to  the  Q  branch  do  not  overlap  each  other  exactly. 

There  is  a  progressive  shift  to  lower  energies  caused  by  the 

vibrational- rotational  interaction.  The  fine  structure  on  the 

(4) 

band  has  been  qualitatively  assigned  by  Lapp  et.  al.  to  be 
due  to  a  large  population  of  an  excited  vibrational  state, 

V  =  1,  whose  transition  is  observed  at  lower  energy  due 

to  anhariHonic  effects .  It  is  obvious  from  these  spectra  that 
the  temperature  dependence  of  the  water  vapor  band  will  not 

be  as  large  as  that  of  N2  and  H2  which  have  been  previously 
investigated  in  flames  .  The  ultimate  sensitivity 
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of  calculated  band  contours  in  comparison  to  experimentally 
obtained  data  is  undoubtedly  limited  by  the  quality  of  the 
experimental  data. 

All  areas  or  intensities  obtained  both  for  cross  section 
calculations  and  temperature  measurements  were  determined 
by  use  of  a  planimeter.  They  represent,  in  most  cases,  an 
average  value  of  ten  determinations  on  each  spectrum. 
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SECTION  III 

experiment/^  data 


3 . 0  GENERAL 

Presented  in  this  section  are  the  high  temperature 
experimental  Raman  spectra  of  H,0  vapor.  All  the  spectra 
were  obtained  using  the  4880  A  laser  line  except  where  cross 

O 

section  scaling  was  carried  out  using  the  5145  A  line.  A 
discussion  of  the  experimental  error  is  also  included  in  this 
section. 

3 . 1  SPECTRA 

The  experimental  apparatus  and  procedure  have  been 
described  in  Section  2.5.  In  Figure  14  are  shown  traces  of 
spectra  of  water  vapor  and  nitrogen  at  various  temperatures 
demonstrating  their  temperature  dependence.  The  spectra  were 
obtained  by  locating  in  the  flames,  with  the  aid  of  a  thermo¬ 
couple,  fairly  homogeneous  temperature  regions.  The  slit  height 
was  adjusted  so  that  only  Raman  scattered  light  from  that  region 
was  collected. 

Included  in  Table  3  are  the  calculated  Relative  Raman 
cross  sections  of  water  vapor  relative  to  nitrogen  for  90® 
scattering.  The  individual  spectra  of  H2O  and  N2  at  the 
various  temperatures,  as  described  in  Table  3,  are  shown 
in  Figures  15  through  18  respectively. 

3.2  SPECTRA  OF  WATER  VAPOR  BETWEEN  1500  K  AND  1750  K 

The  spectra  of  water  vapor  and  N2  at  approximately  1650  K 
and  1750  K  are  shown  in  Figure  19.  The  temperatures  were 
determined  from  the  N2  band  profiles.  Thermocouple  data  at 
these  temperatures  is  unreliable  due  to  an  excess  of  radiation 
cooling  and  catalytic  effects.  Above  1500  K,  the  burners  which 
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FIGURE  18 
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FIGURE  19 

THE  RAMAN  SPECTRA  OF  O  AND  AT  a,  16  50  K.  AND 
1750  K.  THE  SPECTRA  WERE  OBTAINED  WITH  LOCKED-IN 
AMPLIFIER  AND  SYNCHRONOUS  DETECTION  ELECTRONICS. 

THE  EXCITATION  FREQUENCY  WAS  4880  A  WITH  POWER 
LEVELS  OF  NEARLY  8  WATTS.  THE  SPECTRAL  SLIT  WIDTH 

IS  %  8  CM“^ 
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Calculated* 

Average  Thermocouple 

Flow 

Rates ( 1/m) 

90° 

Temperature 

Temperature  (Range) 

“2 

°2 

’^2 

Slit  Heights 

Relative-  Raman 
Cross  Sections 

HjO/Nj 

1100  K 

1107 

(1080-^1145) 

2.150 

1.510 

4.675 

10  mm 

2.47 

1255  K 

1241 

(1239->-1252) 

1.0 

0.5 

2.3 

11.  4 

2.2 

1381  K 

1413 

(1376-'-1440) 

0.5 

0.25 

0.3 

12.7 

2.5 

1488  K 

1500 

(1408->1605) 

0.5 

0.25 

2.0 

8.9 

2.9 

*  Temperature  estimated  from  calculated  N2  spectra. 
+  Corrected  for  radiation  cooling. 
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we  employed  caused  significant  Na  emission  to  result  in  the 
flame  gases,  thereby  making  it  impossible  to  use  the  photon 
counting  detection  electronics.  To  eliminate  the  flame  emission 
signals  from  the  Raman  signals,  we  utilized  a  technique  which 
is  often  used  to  remove  infrared  emission  interferences  from 
infrared  absorption  spectral  data.  The  technique  involves 
chopping  the  source,  in  our  case  the  laser  beam,  and  using 
synchronous  detection  of  the  Raman  signal.  The  method,  in 
effect,  subtracts  the  emission  spectral  data  from  the  observed 
Raman  plus  emission  signals.  The  method  does  have  one  disadvantage 
in  that  considerably  more  noise  is  added  to  the  final  spectrxam. 

With  our  present  experimental  apparatus,  however,  this  was  the 
only  way  we  could  obtain  data.  Table  4  gives  the  temperatures 
and  flow  rates  of  the  pre- flame  gases  and  the  calculated  relative 
Raman  cross  sections  for  this  data. 


TABLE  4.  RELATIVE  RAMAN  CROSS  SECTIONS  OF  H2O  VAPOR 


Temperature  (K) 

Flow  Rates 

(A/m) 

Raman  Cross  Section 
Relative  to  N2  For 
90®  Scattering 
of  4880  A 

*^2 

«2 

°2 

1650 

7.3 

3.  3 

1.3 

2.7 

1750 

7.3 

3.4 

1.7 

■  2.6 

3.3  POLARIZED  RAMAN  SPECTRA  OF  WATER  VAPOR 

The  polarized  and  depolarized  Raman  spectra  of  H2O  vapor  at 
approximately  1075  K  were  obtained  (Figure  20)  by  placing  a 

polarizer  between  the  monochromator  and  the  sampling  area  and 
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FREQUENCY  CM 

FIGURE  20 

THE  POLARIZED  (UPPER  CURVE)  AND  DEPOLARIZED  (LOWER 
CURVE)  RAMAN  SPECTRA  OF  WATER  VAPOR  AT  'v  1075  K. 
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analyzing  the  Raman  scattered  light.  A  polarization  scrambler 
was  placed  immediately  in  front  of  the  entrance  slit  of  the 
monochromator  to  minimize  the  effect  of  the  instruments  selective 
sensitivity  to  polarization.  The  spectra  demonstrated  that  the 
Raman  scattering  of  water  vapor  is  essentially  polarized  and 
results  primarily  from  the  trace  of  the  polarizability  tensor. 

This  greatly  simplifies  the  theoretical  calculations  as 
this  means  that  only  the  Q  branch  of  the  rotation-vibration  mani¬ 
fold  need  be  considered.  This  will  be. discussed  in  more  detail 
in  section  4.6  of  this  report. 

3.4  RAMAN  CROSS  SECTION  OF  H2O  AT  DIFFERENT  WAVELENGTHS 

Listed  below  are  temperature  Relative  Raman  cross  sections 

O  O 

of  H2-O2-N2  flames  using  the  5145  A  and  4880  A  exciting  lines. 

The  flame  temperature  gradients  were  mapped  by  use  of  the  ther¬ 
mocouple  as  in  obtaining  the  previous  data.  The  same  set  up 

O  . 

and  procedure  were  used  as  in  obtaining  data  at  4880  A.  An 
estimate  of  the  temperature  was  also  obtained  using  the  ratio 
of  the  N2  areas.  There  appear  to  be  no  differences  within 
experimental  errors  in  the  Raman  cross  sections  of  H2O  between 
5145  A  and  4880  A. 


Relative 

Raman  Cross  Section  | 

Temperature 

4880  A 

5145  A 

1107  K 

2.47 

2.55 

1488  K 

2.90 

2.70 

3.5  DISCUSSION  OF  EXPERIMENTAL  ERRORS 

There  are  no  significant  errors  in  determination  of  areas 
under  bands  introduced  by  chart  paper  linearity,  monochromator 

I 

and  chart  scan  speed.  Errors  which  arise  from  the  response  of 
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the  chart  pen  to  spectral  events  were  minimized  by  scanning  the 
monochromator  at  rates  much  slower  than  those  of  the  response  of  the 
system  monochromator  and  chart  recorder.  The  uncertainty  in  the 
determination  of  areas  under  the  Raman  profiles  arises  because 
of  the  use  of  a  planimeter.  The  areas  and  their  relative  standard 
deviation  for  10  typical  measurements  are : 

Ajj  Q  =  2565  i  5.94  cps.cm  ^ 

=  2255  -  48.5  cps.cm  ^ 

^2 

In  the  present  study,  a  maximum  of  2%  uncertainty  was  introduced 
into  the  intensity  determination  by  measuring  actual  areas. 

The  errors  in  the  number  densities  of  the  species  of  the 
flame  combustion  products  are  introduced  by  the  flow  rate  and 
temperature  inhomogeneities,  and  determination  techniques. 

The  imprecision  of  the  manufacturer's  calibration  of  the  flow  meters 
is  2%,  or  1  division,  whichever  is  greater.  This  calibration  was 
confirmed  by  our  experimental  calibration.  An  example  of  the  un¬ 
certainties  introduced  by  the  flow  rate  determinations  follows; 

The  flow  rates  of  the  preflame  gases  for  the  1082  K  flame 

were: 

H2  -  2.150  l/m  ±  2% 

N2  --4.675  l/m  ±  2% 

O2  -  1.510  i/m  ±  2% 

Since  for  the  reaction  H2  +  1/2  ^  excess,  the 

amount  of  H2O  produced  is  2.150  i/m  ±  2%  and  the  combustion 
products  of  the  flame  are; 

H2O  -  2.150  l/m  ±  2% 

O2  -  0.435  i/m  ±  2% 

N2  4.675  i/m  ±  2% 
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The  error  in  the  total  number  of  moles  present  in  the  burnt 
flame  can  be  estimated  by  taking  the  square  root  of  the  sum 
of  the  .squares  of  the  individual  gas  uncertainties 
or 

S  =  /(.043)^+  (.0935)^+  (.0302)^ 


S  =  0.10  i/m 


The  relative  deviation,  then,  is  the  uncertainty  divided  by 
the  total  number  of  moles  present,  or 

Relative  Deviation  =  2gQ  =  .014  or  1.4% 


The;  error  in  obtaining  the  mole  fractions  of  the  reaction 
products  (XiT  Xm  )  used  in  calculating  the  Raman  cross 

section  is  obtained  from  their  relative  deviations  and  from 
the  deviations  of  the  total  number  of  moles  present. 


Rel.  Dev. 
Rel.  Dev. 


(H2O) 

(N2) 


2.  5% 
2.5% 


Due  to  the  signal  to  noise  statistics  of  the  spectrum,  there 
will  also  be  systematic  errors  introduced  in  determining  the 
areas  of  the  Raman  bands  1  In  the  worst  case  the  S/N  ratio  is 
20.  This  introduces  an  error  of  at  least  5%  of  the  total  area. 
The  maximum  estimated  error  for  determining  area  ratios  then 
arises  from  the  total  estimated  error  in  each  area  that  is 
to  be  ratioed. 
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Thus,  the  maximum  error  introduced  into  the  relative  density 
determinations  of  H2O  to  NT^  from  flow  rate  and  area  data 
is  : 

Rel.Dev.  /2(5%)^  +  2(2.5%)^  =  8%. 


This  error  estimation  does  not  include  instability  effects  or 
inhomogeneities  within  the  flame. 

The  determination  of  temperature  from  a  ratio  of  the  peak 
intensities  of  the  N2  profile  introduces  not  only  this  error 
but  also  errors  due  to  the  deconvolution  of  the  N2  peaks.  The 
error  in  the  determination  of  the  temperature  from  the  Raman 
nitrogen  spectra  is  a  minimum  of  10%  which  does  not  include 
gradient  effects. 


The  Raman  cross  section  is  given  by 


Raman 


The  standard  relative  deviation  for  „  estimated  from  the 

Raman 

errors  due  to  all  sources  described  above  is 


Rel.Dev.  .  .  =  /(8)^  +  2(2.5)^  =  9% 

Raman ^ 

The  errors  in  the  estimated  temperatures  and  nurdser  densities 
which  arise  from  the  inhomogeneities  in  the  flame  gases ,  such  as 
density  and  temperature  gradients,  are  extremely  difficult  to 
estimate.  Use  of  a  thermocouple  to  find  gases  in  the  flame  with 
the  smallest  temperature  gradients  mimimized  the  temperature 
errors,  and  use  of  burner  designs  which  produced  excellent  flame 
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sovirces  minimized  the  density  fluctuation  errors.  Since 
the  temperature  measurements  and  cross  section  data  are 
obtained  from  ratios  of  band  areas,  these  errors  are  further 
minimized.  They  are  not  completely  cancelled,  however,  since 
the  effects  of  temperature  and  density  gradients,  which  pri¬ 
marily  smear  out  the  band,  are  not  necessarily  the  same  on 
different  Raman  bands. 
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SECTION  IV 

THEORETICAL  CONSIDERATIONS 

4.0  CALCULATION  OF  THE  ROTATIONAL/VIBRATION  SPECTRA  OF  WATER  VAPOR 

4.1  PURE  VIBRATIONAL  SPECTRUM 

4.1.1  The  Harmonic  Oscillator  -  Normal  Coordinate  Vibrational 

Model  ' 

The  equilibrium  geometry  of  water  is  shown  in  Figure  21. 

Water,  a  non-linear  symmetrical  XY2  molecule,  has  symmetry 

Its  elements  of  symmetry  are  C2  axis  through  the  0  atom  and 

bisecting  a  line  between  the  H  atoms  (z-axis) ,  a  plane  of 

symmetry  a  is  determined  by  the  three  nuclei,  and  a  plane  of 
xz 

symmetry  containing  the  C2  axis  is  perpendicular  to  the 

a  .  Non-linear  XY.,  molecules  have  two  normal  vibrations  of 
xz  2 

species  and  V2)  and  one  of  species  and 

are  the  symmetric  and  anti-symmetric  OH  stretches  and  V2  is 
the  bending  mode. 

From  high  resolution  infrared,  far  infrared  and  microwave 
spectroscopy,  the  constants  of  water  molecules  have  been  accu¬ 
rately  determined.  The  fundamental  vibrations  are 

=  3651.7  cm"^ 

a 2  =  1595.0  cm"^  (4-1) 

=  3755. 8  cm  ^ 

The  potential  energy  of  a  harmonic  XY2  molecule  is 

2V  =  +  Arj^)  +  (4.2, 
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2cC 


FIGURE  21 

THE  GEOMETRY  OF  THE  WATER  MOLECULE 
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where 

Ar^^,  Ar2  are  the  changes  of  X-Y  distances, 

6  is  the  change  in  the  apex  angle  2a  (Figure  21) , 
and  are  the  valence  force  constants . 

The  equations  for  the  normal  frequencies  in  terms  of  the 
geometry  and  force  constants  are  as  follows : 

=  (Uy  +  2Ujj  sin^a)K^  (4-3) 

XjL  +  X2  =  4tt2(v^2  +  ^2^) 

=  (Uy  +  Ujj  cos^a)K^  + 

2(Uy  +  2Ujj  sin^a)Kg/il^  (4-4) 

1^X2  =  16ir^Vj^^V2^  =  2Uy(Uy  +  2U^)K^Kg^j^2  (4-5) 

where  is  the  reciprocal  atomic  mass. 

From  equations  (4-3) ,  (4-4) ,  and  (4-5)  and  the  observed 
fundamentals , 

O 

=  7.76  mdyne/A 
=  0.69  mdyne/A 

are  determined. 

The  vibrational  motions  shown  in  Figure  22  are  the  normal 
coordinates -Q^  in  terms  of  which  the- potential  energy,  equation 
(4-2) ,  becomes 

3 

2V  =  ^  X^Q^^  (4-6) 

i=l 
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FIGURE  22 

THE  NORMAL  VIBRATIONS  OF  THE  NON-LINEAR  XY2 


MOLECULE 
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This  approximation  in  polyatomic  molecules  involves  two  separate 
assumptions.  First,  the  complicated  Lissajous '  motion  of  the 
vibrating  nuclei  can  be  divided  into  simple  motions  where  each 
nucleus  is  moving  in  a  straight  line  through  its  equilibrium 
position.  And  second,  these  motions  are  harmonic  oscillations. 
The  normal  coordinates  shown  in  equation  (4-6)  are  related  to 
the  internal  coordinates,  q^,  shown  in  equation  (4-2)  by  a 
linear  transformation,  i.e. 

3 

■’1=  L  ‘ikOk 

K=1 

Only  the  internal  force  constants  K^,  shown  in  equations  (4-3)  , 
(4-4) ,  and  (4-5) ,  have  real  meaning  since  they  depend  on  the 
local  bonding.  They  can  be  transferred  from  one  molecule  to 
another  where  identical  or  even  similar  bonding  situations 
occur,  and  they  are  the  basis  for  characteristic  vibrational 
frequencies . 

For  a  harmonic  oscillator,  the  energy  values  are  given 
by 

Ei  =  hv^(n^  +  1/2),  n^  =  0,  1,  2,  -  (4-8) 

where 


The  total  vibrational  energy  for  triatomic  molecules  is 

G(nj^,n2n2)  =  hf  =  +  1/2)  +  02(n2  +  1/2) 

+  a3(n3  +  1/2)  (4-9) 
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where 


"i  =  '’i/c 


are  the  classical  vibrational  frequencies  measured  in  cm 
units.  The  vibrational  quantum  numbers  are  n^^,  n2f 


4.1.2  The  Anharmonic  Vibrational  Model 


When  the  cubic  and  quatric  terms  are  considered,  the 
potential  energy  equation  in  normal  coordinates  becomes 

3  3  3  3 


2V 


i  i  j  K 


2^  1^2^ 2^ 

i  j  K  A 


(4-10) 


Equation  (4-9) ,  the  vibration  energy  equation  for  a  non-linear 
triatomic  molecule,  becomes 


G(nj^,n2,n2)  =  U]^(n^  +  1/2)  +  W2  (1^2 

+  0)3  (n3  +  1/2)  +  X^j^(n^  +  1/2)^ 

+  X22tn2  +  1/2)^+  X33(n3  +  1/2)^ 

1 

+  Xj^2 

+  Xj^3(n3^  +  1/2)  (n3  +  1/2) 

+  X23  (n2  +  1/2) (n3  +  1/2)  (4-11) 
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The  are  the  frequencies  for  Infinitesimal  amplitudes  or 

zero-order  frequencies  and  were  first  approximated  by  the 

in  equation  (4-9).  The  are  the  anharmonicity  constants. 

The  X . .  like  the  u  X  used  for  diatomic  molecules  are  a  measure 
11  e  e 

of  the  anharmonicity  but  the  x^j^(i^K)  are  a  measure  of  the 
interaction  between  the  vibrations.  When  the  fundamental 
is  excited  the  energy  of  vibration  is 

=  G(1,0,0)  -  G(0,0,0) 

=  +  2X^^  +  1/2X^2  V2Xj^3  (4-12) 


An  interpretation  of  equation  (4-12)  is  that  when  is 
excited  most  of  the  energy  stimulates  an  enharmonic  motion 
but  a  small  amount  of  it  also  stimulates  the  motions  Q2 
and  Q3  (see  Figure  22)  . 

The  constants  for  H2O  are  (all  in  cm”^) : 


=  3825.32  =  -43.89  Xj^2  ==  “20.02 

0)2  =  1653.91  X22  =  -19.5  X^3  =  -155.06 

0)3  =  3935.59  X33  =  -46.37  X23  =  -19.81 

(4-13) 

There  is  a  very  large  interaction  (155  cm  ^)  between  the 
symmetric  and  antisymmetric  stretching  vibrations.  Comparing 
cr.  (observed  band  centers)  with  (the  zero-order  frequencies) 
a  measure  of  the  anharmonicity  of  the  vibrations  is  derived. 


(4-14) 


56 


AEDC-TR-75-1 50 


then 

i  !!i 

1  4.5% 

2  3.6% 

3  4.6%  (4-15) 

With  equation  (4-11)  and  the  constants  shown  in  equation 
(4-13) ,  we  can  use  energy  levels  for  all  excited  states  and, 
by  differences  between  them,  calculate  the  observcdole  vibrational 
transitions.  The  selection  rule  for  both  infrared  and  Raman 
transitions  in  non-linear  XY2  molecules  is  An^  =  +1.  This  is  a 
harmonic  oscillator  approximation  rule,  but  since  the -vibra¬ 
tions  are  more  than  95%  harmonic  it  is  relatively  accurate.  At 
temperatures  well  above  room  temperature,  the  equilibrium  popu¬ 
lation  of  excited  states  (energy  levels)  will  be  significant. 
Vibrational  bands  in  which  the  initial  state  is  other  than  the 
ground  vibration  state  are  known  as  difference  bands  whose 
frequency  will  vary  from  the  fundamental  by  the  relative 
anharmonicity  of  the  involved  states.  These  bands  will  have 
the  general  formula 

t^h  “  ^2'  ”3^  “  ^2*  (4-16) 

Since  the  anharmonicity  is  small,  the  transition  probability 
of  these  bands  will  be  identical  to  that  of  the  fundamental. 

Thus,  the  intensity  ratio  of  a  difference  band  to  the  funda¬ 
mental  is  very  nearly  equal  to  Boltzmann  factor,  expl- (hco^/kT) ] , 
which  gives  the  relative  population  of  the  states  at  a  given 
temperature.  The  relative  intensities  of  the  hot  band  to  the 
fundamental  for  the  region  are  shown  in  Table  5.  In  the 
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Raman  spectrum  of  water  vapor,  only  the  fundamental  should 
be  seen  at  temperatures  below  1000  K;  from  1000  K  to  2000  K  a 
second  line  centered  20  cm  ^  red  shifted  from  the  fundamental 
should  be  clearly  visible.  Above  2000  K,  three  more  bands 
should  become  observable. 

4.1.3  Resonant  Interactions 

When  two  vibrational  levels ,  only  one  of  which  can  be  a 
fundamental,  have  the  same  symmetry  and  have  nearly  the  same 
energy,  that  is,  may  be  accidentally  degenerate,  they  can  inter¬ 
act  to  produce  a  perturbation  called  Fermi  resonance.  The 
result  of  such  an  interaction  is  that  the  levels  "repel”  each 
other  and  there  is  a  mixing  of  intensity.  In  water  this  occurs 
between  states  having  the  following  quantum  nvimber  relations: 

{n^  +  2,  TL^,  n^)  and  (n^^,  +  '2).  Of  the  states  shown 

in  the  table  only  (200)  and  (210)  are  effected  states.  In 
these  cases  the  perturbations  are  small  and  these  states  do 
not  contribute  significantly  to  the  observable  spectrum. 

4.1.4  Fermi  Resonance 

In  the  previous  section,  the  strong  resonant  interaction 
between  2v^(7201.40  cm  ^)  and  2^2^7445.00  cm  ^)  was  shown  not 
to  contribute  to  the  observed  spectrum.  There  is  also  a  weak 
interaction  between  (3657. 054)  and  2v2 (3151. 631)  .  For  Fermi 

resonance  between  two  energy  levels , 

=  4/W^2/^  (4-17) 

where  A  is  the  perturbed  observed  separation,  5  is  the  unper¬ 
turbed  separation,  and  W^2  is  the  energy  of  interaction  and 
is  related  to  the  enharmonic  force  constant  ^2.22 
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VIBRATIONAL  TRANSITIONS  AND  RELATIVE  INTENSITIES 
AT  VARIOUS  TEMPERATURES  IN  THE  REGION 


1  Transition 

Lower 

State 

Energy 

(cm”-*-) 

Upper 

State 

Energy 

(cm”^) 

Band 

Center 

(c.m”-*-) 

Relative  Intensity  * 

Lower 

State 

Quantum 

Numbers 

Upper 

State 

Quantum 

Numbers 

300 

(K) 

1000 

(K) 

1250 

(K) 

1500 

(K) 

1750 

(K) 

2000 

(K) 

3000  i 

(.K) 

(000) 

(100) 

0* 

3650.0 

3650.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000  1 

(010) 

(110) 

1590.0 

5225.0 

3630.0 

0.000 

0.101 

0.160 

0.217 

0.269 

0.317 

0.466  * 

(020) 

(120) 

3151.0 

6761.0 

3610.0 

0.000 

0.011 

0.027 

0.049 

0.075 

0.10-1 

0.221  j 

!  (100) 

(200 

3650.0 

7212.2 

3562.2 

0.000 

0.005 

0.035 

0.030 

0.050 

0.072 

0.1/4 

(001) 

(101) 

3755.4 

7250.4 

3495.0 

0.000 

0.005 

0.013 

0.027 

0.046 

0.067 

0.165  i 

(030) 

(130) 

4668.0 

8257.9 

3589.9 

0.000 

0.001 

0.005 

0.011 

0.022 

0.032 

0.101  1 

(110) 

(210) 

5225.0 

8767.2 

3542.2 

0.000 

0.001 

0.002 

0.007 

0.014 

0.023 

O.CSl  ' 

(011) 

(111) 

5330.6 

8305.5 

3474.9 

0.000 

0.000 

0.002 

0.006 

0.012 

0.022 

Q.07S 

(040) 

(140) 

6146.0 

9715.9 

3569.9 

0.000 

0.000 

0.001 

0.003 

0.006 

0.012 

0.C52  ; 

*  Assigned  value 
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TABLE  6 


GROUND  VIBRATIONAL 

ROTATIONAL  ENERGY  LEVELS 

STATE 

FOR  J=0  TO 

5 

E  (cra”^) 

Relative  Population  at 

300  K  1000  K  200C  K 

0(0,0) 

0 

1.00 

1.00 

1.00 

1(0,1) 

23.8 

0.89 

0.97 

0.98 

1(1,1) 

37.1 

0.84 

0.95 

0.97 

1(1,0) 

42.3 

0.82 

0.94 

0.97 

2(0,2) 

70.2 

0.71 

0.90 

0.95 

2(1,2) 

79.5 

0.68 

0.88 

0.94 

2(1,1) 

95.2 

0.63 

0.87 

0.93 

2(2,1) 

135.0 

0.52 

0.82 

0.91 

2(2,0) 

136.0 

0.52 

0.82 

0.91 

3(0,3) 

136.9 

0.52 

0.82 

0.91 

3(1,3) 

142.9 

0.52 

0.82 

0.91 

3(1,2) 

170.6 

0.44 

0.78 

0.90 

3(2,2) 

206.4 

0.37 

0.74 

0.86 

3(2,1) 

232.7 

0.33 

0.72 

0.85 

3(3,1) 

286.6 

0.25 

0.66 

0.81 

3(3,0) 

286.8 

0.25 

0.66 

0.81 

4(0,4) 

222.3 

0.34 

0.73) 

0.85 

4(1,4) 

225.0 

0.34 

0.73 

0.85 

4(1,3) 

275.9 

0.27 

0.67 

0.82 

4(2,3) 

300.6 

0.24 

0.65 

0.81 

4(2,2) 

316.9 

0.22 

0.63 

0.80 

4(3,2) 

383.8 

0.16 

0.58 

0.76 

4(3,1) 

385.2 

0.16 

0.58 

0.76 

4(4,1) 

484.6 

0.10 

0.50 

0.71 

4(4,0) 

487.7 

0.10 

0.50 

0.71 

5(0,5) 

325.9 

0.21 

0.63 

0.79 

5(1,5) 

327.1 

0.21 

0.63 

0.79 

5(1,4) 

400.4 

0.15 

0.56 

0.75 

5(2,4) 

416.9 

0.14 

0.55 

0.74 

5(2,3) 

447.8 

0.12 

0.53 

0.72 

5(3,3) 

505.6 

0.09 

0.48 

0.69 

5(3,2) 

510.7 

0.09 

0.47 

0.69 

5(4,2) 

615.2 

0.05 

0.41 

0.64 

5(4,1) 

615.5 

0.05 

0.41 

0.64 

5(5,1) 

758.5 

0.03 

0.34 

0.58 

5(5,0) 

758.5 

0.03 

0.34 

0.58 
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TABLE  7 

THE  CONSTANTS  FOR  THE  WATER  MOLECULE  (cm~^) 


Vibration 

1 

B 
a  • 

1 

C 

1 

+0.495 

+0.224 

+0.145 

2 

-2.659 

-0.202 

+0.105 

3 

+1.234 

+0.112 

+0.169 

TABLE  8 

ROTATIONAL  CONSTANTS  FOR  THE  GROUND 
AND  GENERAL  EXCITED  STATES  OF  WATER  (cm~^) 


Vibrational 

State 

A 

B 

C 

b 

000 

27.795 

14.508 

9.289 

-0.1642 

100 

27.300 

14.284 

9.144 

-0.1649 

010 

30.454 

14.710 

9.184 

-0.1492 

110 

29.959 

14.486 

9.039 

-0.1497 

020 

33.113 

14.912 

9.079 

-0.1381 

120 

32.618 

14.688 

8.934 

-0.1383 
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”1^  ^  <^2"^  2)^/^  (4-18) 

Since  k^22  been  found  to  be  71.71  cm  then  bhe 

interaction  of  with  2^2  is  35.86  cm  Then  for  a  A  = 

505.423  cm~^,  5  =  500.310  cm”^.  The  Fermi  shift  (F  )  is 

a 

F  =  (A-5)/2  =  2.557  cm"^  (4-19) 

3i 

The  unperturbed  frequency  (v^^)  ®  is  3654.  497  cm 

Since  the  observed  band  centers  are  already  known  from 
high  resolution  infrared  work,  the  Fermi  resonance  calculation 
is  really  incidental  to  the  problem,  which  is  to  calculate 
the  observed  rotation- vibration  band  envelope  in  the  region. 

As  shown  in  the  previous  section,  the  only  hot  band  signi¬ 
ficantly  contributing  to  our  spectrum  is  (v^+V2“V2^  (010)  ■* 

(110).  Now  (110)  and  (010)  have  been  observed  at  5235.00  and 
1594.59  cm  ^  respectively.  The  observed  (perturbed)  hot  band 
taking  the  difference  is  at  3640.41  cm~^.  The  unperturbed 

band  center  of  the  hot  band  with  W.-=  62.10  and  A  =  568.30  cm  ^ 

-1 

rs  at  3626.67  cm  .  Due  to  the  larger  interaction  constant,  the 
Fermi  shift  is  13.74  cm  ^  in  the  hot  band  versus  2,56  cm”^ 

in  the  fundamental.  For  our  purposes,  use  of  =  3657.05  and 

Vj^+V2-V2  =  3640.41,  the  observed  band  centers  inherently 

includes  corrections  for  Fermi  resonance.  The  frequencies 

presented  here  vary  slightly  from  those  used  initially  since 

we  are  now  using  a  set  of  frequencies  and  constants  compiled 

(5) 

by  D.  Smith  and  J.  Overend '  ,  which  we  believe  are  the  most 

up-to-date  and  accurate  of  those  published.- 
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4.1.5  The  Rotational  Structure  of  Vibrational  Bands 

Our  discussion  so  far  pertains  to  the  calculation  of 
the  vibrational  band  centers.  About  each  band  center  is  a 
rotational  manifold  representing  transitions  having  the  same 
change  in  vibrational  quantum  numbers  but  different  in  their 
rotational  quantum  numbers. 

4. 1.5.1  Asymmetric  Top  Rigid  Rotor 

Using  the  geometry  shown  in  Figure  21  and  assigning 
the  molecule  fixed  axis  XYZ  according  to  Section  4.1.1,  the 
instantaneous  moments  of  inertia, 

N 

^XX  "  ,  etc.  (4-20) 

i=l 

are  =  1.02024  x  lO"*^®  g  cm^,  1^^  =  2.93786  x  10"^°,  and 

-40 

^ZZ  “  ^  .91762  X  10  .  The  instantaneous  products  of  inertia 

N 

IjjY  ~  ^i^i^i  ^  etc.  (4—21) 

i=l 

are  1^^  =  =  0.  When  the  cross  terms  of  the  inertial 

tensor  are  zero,  the  diagonal  terms  are  the  principal  moments 
of  inertia.  When  all  three  moments  of  inertia  are  different, 
the  molecule  is  called  an  asymmetric  top.  It  is  customary  in 
molecular  spectroscopy  to  designate  these  axes  as  a,  b,  c 
with  Iq*  It  is  also  convenient  to  use  the  rotational 

constants 

A  =  hV2cI^,  etc.  (4-22) 
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with  A  >  B  >  C.  For  water  the  rotational  constants  are  (all 
in  cm  : 


A  =  27.4355 
e 

II 

27.795 

B  =  14.5966 
e 

II 

o 

CQ 

14.508 

(4-23) 

C  =  9.5276 

e 

^O  = 

9.289 

The  first  set  of  constants  in  equation  (4-23)  are  for 
the  equilibrium  configuration  while  the  second  set  are  for 
the  ground  vibrational  state.  Their  difference  will  be  dis¬ 
cussed  later.  The  energy  matrix  for  the  rotational  states  with 
the  quantum  numbers  J  and  K  and  using  the  following  notation 
for  matrix  elements: 

(J,k|Hj^|  J' ,K')  =  y*i|;*(J,K)Hj^ijjj^(J' ,K')dT  (4-24) 

where  for  the  rigid  rotor  the  Hamiltonian  is 


The  matrices  are 

(J,K|Hj^|  J,K)  =  1/2(Aq  +  C^)IJ(J  +  1)  -  K^]  +  2BqK^ 

.  (J,k|Hj^|J,K  ±  2)  =  1/4(Cq  -  A^)  |[J(J  +  1)  -  K(K  -  1)  ] 

[J(J  +1)  -  (K  ±  1)  (K  i  2)  ] 

(4-26) 


where  K=-J,  ...,  0,  ...»  J 
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Equation  (4-26)  will  give  a  symmetric  matrix  of  the  following 
type  for  each  value  of  J. 


X| 

-J 

-J+1 

-J+2 

-J+3 

J 

-J 

^0 

0 

0  _ 

•  •  • 

-.T+1 

0 

^1 

0 

Id  2  •  •  • 

«  •  • 

-J+2 

*>1 

0 

^2 

0 

•  •  • 

-J+3 

0 

CM 

Si 

0 

K3  ... 

•  •  • 

J 

•  *  • 

•  •  • 

•  •  • 

•  •  ■  »  •  • 

«J+1 

(4-27) 


To  simplify  calculations  this  matrix  can  always  be  trans¬ 
formed  (block  diagonalized)  into  four  submatrices  labeled 

E  ,  0^,  O”  where  E^  means  having  an  element  where  K  is  even 
and  positive.  Diagonalization  of  these  matrices  gives  the 
rotational  energy  levels.  A  program  was  developed  to 
calculate  these  energy  levels  for  each  J. 

The  Wang  asymmetry  parameter,  (often  used  to  describe 
vibration-rotation  spectra) 

^  ~  2  I A- 1/2 (B+Cj ]  (4-28) 

is  a  measure  of  a  molecule's  asymmetry.  For  the  equilibrium 
constants  in  equation  (4-23),  b  =  0.1649.  In  the  prolate  limit, 
A>B=C,b=0  and  in  the  oblate  limit,  A=B>C,  b=-l. 
Typical  asymmetric  energy  levels  for  low  J  are  shown  in  Figure 
23.  Note  that  it  is  best  to  show  water  as  a  perturbation 
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PROLATE  TOP 


HoO 


OBLATE  TOP 


k 

-IjO 

-0.5 

0.0 

0.5 

1.0 

b 

0.0 

-0.142 

-0.333 

-0.6 

-1.0 

FIGURE  23 

ASYMMETRIC  ENERGY  LEVELS  FOR  LOW  J 
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of  the  prolate  limit  rather  than  of  the  oblate  limit.  Another 
way  to  designate  the  energy  levels  is  the  bipartite  index 
K_iKi  where  these  are  the  energy  levels  in  the  prolate  and 
oblate  limits  respectively.  Also 

K  =  K_^  -K^  (4-29) 

The  energy  levels  for  J  =  0  to  5  are  shown  in  Table  6  along 
with  the  relative  population  at  various  temperatures.  Unlike 
those  intensity  values  in  Table  5,  the  transition  proba¬ 
bilities  for  rotation  are  not  constant  over  J  and  K_j^.  There¬ 
fore,  the  relative  population  of  the  levels  is  not  a  measure 
of  the  relative  intensity.  The  rotational  line  strengths  will 
be  considered  later;  however,  they  increase  more  rapidly  with 
K_^  than  with  J.  The  low  levels  of  high  J  contribute 
extensively  to  the  observed  band  shape. 

4. '2  CALCULATION  OF  THE  ROTATION-VIBRATION  SPECTRA  IN  THE 
REGION  OF  WATER  VAPOR 

4.2.1  Higher  Order  Approximation  to  the  Energies  of  an 
Asymmetric  Rotor  -  Artificial  Distortion 

The  effects  due  to  centrifugal  distortion  are  usually  too 
small  to  observe  at  frequency  accuracy  and  resolutions  of  less 
than  0.1  cm  The  treatment  of  the  non-rigid  rotor  adds  terms 
of  the  following  type  to  the  energy  matrix:  (K/K) ,  (K/K±2) ,  and 
(K/Ki4) .  A  better  (should  converge  more  rapidly)  equation  for 
the  diagonal  elements  than  the  one  shown  previously  includes 
the  centrifugal  stretching  constants; 
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E  =  2+2  J{J+1)  +  ^  A  - 

'  '  (4-30) 

-  DjJ^  (J+1)^  -  Djj^  J(J+1)  -  Dj^ 

The  last  terms  in  equation  (4-30)  are  due  to  centrifugal 

distortion.  It  is  easily  observed  that  these  terms  become 

significant  at  high  J  and  K.  The  exact  values  for  H-, 

and  D  have  not  as  yet  been  determined,  but  they  are  on  the 
^  “  4  5  “1 

order  10  to  10  cm  .  But  as  we  have  shown,  the  selection 
rules  are  AJ  =  0  euid  AK  =  0.  Hence,  the  shift  in  spectral 
lines  should  have  the  form: 

(D^  -  Dq)  (J+1)  ^+J(J+1)K^+  (4-31) 

where  D  .  ~  change  in  centrifugal  distortion  between 

^  "■  6  —  1 
the  groimd  state  and  v^,  is  at  most  of  the  order  10  cm  .  For 

example,  a  line  that  should  contribute  to  our  observed  spectrum 

is  J  =  20  and  K  =  9,  is  shifted  only  0.2  cm 

4.2.2  Higher  Approximations  to  the  Energies  of  an  Asymmetric 
Rotor 


4. 2. 2.1  Rotation-Vibration  Interaction 

The  rotational  constants  of  excited  vibrational  states 
will  be  different  from  those  in  the  ground  state  since  the 
geometry  of  the  molecule  changes  with  excitation.  This  depend¬ 
ency  has  been  determined,  using  harmonic  oscillator  wave  fimctions, 
and  is  as  follows: 

A  =  A  -  aL,^(n.  +  1/2)  ,  etc.  (4-32) 

V  e  1  1 
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Since  this  is  a  small  perturbation  to  the  rigid  rotor  energy 
levels,  the  harmonic  approximation  is  adequate.  For  water 
the  constants  a .  are  shown  in  Table  7 .  The  rotational  con¬ 
stants  for  the  vibrational  states  involved  in  the  fundamental 
and  the  two  most  intense  hot  bands  as  calculated  in  Table  5 
are  tabulated  in  Table  8.  The  tabulation  of  the  Wang  asymmetry 
parameter  shows  there  is  very  little  change  in  asymmetry  for 
An^^  =  1,  and  hence  this  does  not  contribute  to  the  band  struc¬ 
ture. 

4.3  CORIOLIS  PERTURBATION 

Two  vibrations  may  interact  through  Coriolis  coupling 
when  the  "product"  of  their  species  contains  the  species  of 
a  rotation  of  the  molecule  (Jahn's  rule).  The  interaction  is 
significant  only  when  the  two  vibrations  are  nearly  degenerate. 
For  water  at  .3652  cm”^  and  at  3756  cm~^  are  close  enough 
to  strongly  interact  and 

X  Bj,(V3)  =  Bj^(Ry).  (4-33) 

Since  Y  has  been  identified  with  the  rotation  constant  C,  this 
will  be  C-type  Coriolis  Coupling  whose  nonvanishing  matrix 
elements  are 

=  i  f(J,K,+  (4-34) 

=  -i  f(J,K,+  )Rj^^  (4-35) 


i  =  /=! 

and 

f(J,K,+  )  =  %[J(J  +  1)  -  K(K  +  1)]^ 


and 


I  V3,J,K+1) 


(V3,J,k1Hj^'  I  Vj^,J,K+1) 


where 


C  Y 

V  =  =13 


*  “3 

(ujj^W3) 
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TABLE  9 

CALCULATED  AND  OBSERVED  Q  BRANCH  TRANSITIONS  BETWEEN 
THE  VIBRATIONAL  STATES  (000)  AND  (100)  FOR  H^O  VAPOR 


State 

Observed 

0(0,  0) 

3657.054 

1(0,  1) 

3656.66 

1(1,  1) 

3656.171 

1(1,  0) 

3656.114 

2(0,  2) 

3655.85 

2(1,  2) 

3655.41 

2(1,  1) 

3655.29 

2(2,  1) 

3653.79 

2(2,  0) 

3653.81 

3(0,  3) 

3654.62 

3(1,  3) 

3654.26 

3(1,  2) 

3655.03 

3(2,  2) 

3652.58 

3(2,  1) 

3652.62 

3(3,  1) 

3649.98 

3(3,  0) 

3649.93 

4(0,  4) 

3652.95 

4(1,  4) 

3652.73 

4(1,  3) 

3652.28 

4(2,  3) 

3650.96 

cm  ) _ Calculated  (cm  ) 


3657.054 

3656.695 

3656.414 

3656.335 

3655.962 

3655.755 

3655.518 

3654.691 

3654.705 

3654.907 

3654.774 

3654.302 

3653.598 

3653.531 

3652.037 

3652.036 

3653.553 

3653.480 

3652.699 

3652.133 
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TABLE  9.  (Continued) 


4(2, 

2) 

3650.77 

3651.943 

4(3, 

2) 

3648.34 

3650.537 

4(3, 

1) 

3648.02 

3650.529 

4(4, 

1) 

3646.86 

3648. 385 

4(4, 

0) 

3646.64 

3648. 385 

5(0, 

5) 

3650.95 

3651.916 

5(1, 

5) 

3650.83 

3652. 880 

5(1, 

4) 

3650  .08 

3650.731 

5(2, 

4) 

3648.94 

3650.318 

5(2, 

3) 

3649.41 

3649.921 

5(3, 

3) 

3646.36 

3648. 658 

5(3, 

2) 

3645.12 

3648. 621 

5(4, 

2) 

3647.65 

3646.509 

5(4, 

1) 

3645.89 

3646.509 

5(5, 

1) 

3640.00 

3643.743 

5(5, 

0) 

3639.81 

3643. 743 
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is  the  Coriolis  coupling  constant.  For  each  J,  there 
will  be  a  Hermetian  matrix  (a^^^  =  af where  *  means  complex 
conjugate)  of  order  4J  +  2.  If  are  the  2J  +  1  matrices 
described  in  Section  4. 1.5.1,  the  block  diagonalized  energy 
matrix  will  have  the  form 


(4-36) 


where  u^2  is  2J  +  1  matrix  with  the  elements  described  in 
equation  (4-34).  This  matrix  can  be  further  block  diag¬ 
onalized  into  four  submatrices  ofJ+1,  J+1,  j  and  J 
dimensions  where  the  interacting  sets  of  rotation-vibration 
states  will  be;  E  *  »  0^  and  E  ♦-►O  . 


Fraley  and  Rao^®^  give  the  observed  energy  levels  of 
H2O  up  to  J  =  9  for  the  vibrational  states  (000) ,  (100) ,  and 
(001) .  We  have  calculated  the  ground  state  levels  for  up  to 
J  =  5  without  including  the  effects  of  Coriolis  interactions. 

The  comparison  of  the  observed  and  calculated  results.  Table  9, 
shows  maximum  differences  of  no  greater  than  1  to  2  cm  Thus, 
we  can  assume  that  Coriolis  coupling  can  be  neglected  for  low 
values  of  J.  For  large  values  of  J,  J  >  15,  this  interaction 
may  yet  have  to  be  included  in  the  calculations.  This  will  be 
determined  when  the  initial  calculations  are  carried  out  for 
temperatures  greater  than  1000  K. 
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4.4  SELECTION  RULES  FOR  A-TYPE  BANDS  IN  AN  ASYMMETRIC  TOP 

Once  the  energy  levels  of  the  upper  and  lower  states  have 
been  determined,  the  selection  rules  tell  us  what  transitions 
can  occur.  The  selection  rules  for  Raman  transitions  where 
only  ttyy,  and  are  different  from  zero  are: 

A J  =  0,  ±1,  ±2 

ee  ee 
eo  eo 
oe  oe 

oo  oo  (4-37) 

where  ee  etc,  refer  to  the  bipartite  notation  means 

both  K_^  and  are  even.  The  dominant  transitions  are  those 
where  AJ  =  0,  the  so  called  Q  branch.  This  sharp  strong 
central  "line"  will  in  general  be  the  only  feature  that  is 
observed.  The  changes  AK_j^  and  aK^^  can  beO,  ±2,  ±4,  .  .  . 
However,  those  with  AK_j^  and  AK^^  =  0  will  be  dominant.  Con¬ 
sequently  the  rotational  structure  about  and  its  hot  bands 
will  be  primarily  transitions  where 

^1'  +  1,  J(K_^,  K^).  (4-38) 

There  are  two  reasons  why  these  lines  do  not  exactly  coincide: 

(a)  the  asymmetry  of  the  molecule  changes  on  excitation 
of  v^.  This  was  shown  in  Table  8  to  be  relatively 
small. 

(b)  Coriolis  coupling  with  v^. 

4.5  CALCULATION  OF  RELATIVE  LINE  STRENGTHS 

The  transition  probability  for  Raman  lines  is 

I  (V,r1  P®|  V  ,R' )'l  ^  (4-39) 
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where  P**  =aE.  P**,  the  time- independent  induced  dipole  moment 
is  equal  to  the  polarizability,  a,  times  the  applied  field. 

IV’fR')  1^  (4-40) 

g  g' 

It  is  possible  to  separate  the  rotation  and  vibration  parts  of 
the  integral  in  equation^ (4-40)  as  follows 


(V,R|agg,lV',R-)  =  (Vlagg,|v')(R|<DpgCt)j.g,lR') 

(4-41) 

It  is  not  necessary  to  know  the  exact  values  of  the  (v|a  ,|v'), 
merely  that  the  quantity  is  nonvanishing.  Hence  the  line 
strengths  can  be  defined  as 


(4-4  2) 


and,  by  combining  these  appropriately  with  the  lower  state 
Boltzmann  population,  the  relative  intensities  of  vibrational 
rotational  transitions  may  be  calculated.  The  Coriolis 
coupling  leads  to  an  intensity  perturbation  also  and  must  be 
taken  into  account.  The  eigenvectors  determined  in  the  diagonal 
ization  of  the  4J  +  2  Hermitian  matrix  can  be  used  to  find 
the  degree  of  character  in  the  line.  Since  in  the  Raman  spec 
trum  the  intensity  of  is  nearly  zero,  multiplication  by  one 
minus  this  factor  will  be  used  to  reduce  the  intensity. 

4.6  CALCULATION  OF  THE  BAND  ENVELOPE 

Once  the  intensity  and  frequency  of  the  component  lines 
have  been  determined,  the  overall  band  envelope  at  a  given 
temperature  can  be  determined. 
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Gaufres 
Raman  line 


gives  for  the  intensity  of  a  vibration-rotation 


pCT 

(V 

,  J', 

n' , 

M' 

^  V”,  J",  n",  M")  a 

|<v"| 

a.  , 

3 

|v- 

>1^  |<  J",  n",  M"|$ 

pa3k 

j» 

k 

(4-43) 


A  vibrational  band  is  allowed  when  at  least  one  of  the  six 
different  elements  of  <V"ta.,  lv'>  is  non-vanishing,  while  the 
rotational  structure  of  an  allowed  band  arises  from  the  non¬ 
vanishing  elements  of  the  second  term  of  the  product.  In  ex¬ 
pression  4-43,  the  molecular  coordinates  are  represented  by  p 
and  a,  the  laboratory  coordinates  by  j  and  k,  the  single  prime 
and  double  prime  denotes  the  upper  and  lower  states  of  the  vi¬ 
brational,  rotation  and  magnetic  quantum  states  respectively 
for  V,  J,  M  and  n  where  n  represents  a  rotational  quantum 
number  other  than  J.  The  term  is  the  polarizability  tensor 
■of  the  molecule  and  the  orientation  matrix  of  the 

molecule  in  the  laboratory  reference  frame.  Its  elements  are 
the  Euler  angles  defined  by  the  system. 

The  polarizability  tensor  may  be  written  as 


a 

XX 
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a 
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a 

yy  yz 
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in  which 


a=-=-(a  +  a  +  a) 

3  XX  yy  z  z 


(4-44) 


and 


and 


“  jj  =  "jj  -  “ 


“'jk  =  “jk  'J 
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The  trace  of  the  LHS  of  equation  4-44  is  given  by  the  first 
term  of  the  RHS  of  equation  4-44  and  the  anisotropy  of  a  is 
given  by  the  second  term  in  the  RHS.  The  selection  rules  for 

(7) 

trace  scattering  are  known  to  be: 


J*  =  J",  n*  =  n",  and  M'  =  M". 

Thus,  trace  scattering  can  only  be  in  the  Q  branch  (AJ  =  0), 
whereas  the  anisotropy  scattering  occurs  from  Or  P,  R,  and  S 
branches  (AJ  =  ±1,  ±2} . 

The  trace  scattering  and  the  anisotropy  scattering  may  be 
separated  experimentally  since  the  polarized  spectra, 
the  dipolar ized  spectra,  are  functions  of  the  trace  and 

anistropy  scattering: 


I 


a  451  +  41 

"  tr  anis 


(4-45) 


“  ”anis 

We  have  shown  experimentally  that  the  water  vapor  spectrum  may 
be  safely  approximated  as  only  trace  scattering. 

Gaufres  gives  for  the  line  intensity  of  trace  scattering 
for  asymmetric  tops  the  following  expression: 


I(J)  =  P(J)  X  (2  J  +  Dg"  exp(-E"^^/kT)  (4-46) 

In  this  expression  the  pre-exponential  term  F ( J)  is  approximated 
to  be  unity.  This  is  only  valid  for  isotropic  scattering?  for 
anisotropic  scattering  F(J)  becomes  a  complicated  function  depen- 
dnet  on  symmetry.  The  term  is  the  ground  state  energy  level 

of  the  transition.  The  theoretical  evaluation  of  it  has  been  dis¬ 
cussed  previously.  The  term  g"  is  the  nuclear  spin  degeneracy  in 

the  initial  state.  For  the  water  molecule  it  has  a  value  of  3 

( 8 ) 

for  symmetric  states  and  1  for  anti-symmetric  states. 
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SECTION  V 

CALCULATED  EAMAN  BAND  PROFILES  OF  THE  REGION  OF  H2O  VAPOR 

5.  1  COMPUTER  PROGRAM  FOR  CALCULATION  OF  THE  Q  BRANCH  OF  H2O 
VAPOR  IN  THE  REGION 

The  problem  of  calculating  the  Raman  band  profile  of  water 
vapor  in  the  region  is  greatly  simplified  by  the  absence  of 
any  significant  contributions  from  anisotropic  scattering.  This 
fact  was  clearly  demonstrated  in  Section  3.3  in  the  discussion 
of  the  polarized  Raman  spectra.  The  selection  rules  for  the 
water  molecule  predict  that  only  the  Q -branch  of  the  rotation- 
vibration  transition  manifold  may  be  active  in  isotropic  Raman 
scattering. 

5.1.1  Limitations  and  Description  of  Calculations 

A  computer  program  was  generated  (Appendix)  to  calculate 

the  Raman  band  profile  of  H2O  vapor.  The  program  was  written  in 

such  a  manner  as  to  design  a  theoretical  model  which  reproduced 

the  experimental  observation  as  closely  as  possible.  For  exaii^>le, 

the  experimental  spectra  were  obtained  with  a  spectral  slit  width 
-1  " 

of  8  cm  (=:3  A)  and  the  conputer  program  uses  a  slit  function 
with  the  same  width.  No  provisions  were  made  for  changing  this 
function  with  external  input  statements. 

Since  the  calculation  involves  analyzing  very  low  resolution 
spectra,  the  perturbation  of  line  position  due  to  centrifugal 
distortion  and  Coriolis  coupling,  and  the  effects  of  pressure 
broadening  doppler  shifts  and  natural  line  widths  were  not 
included  in  the  calculation  of  the  contour  shapes.  These  effects, 
which  perturb  the  line  widths  of  the  Q  branch  of  H2O  by  tenths 
of  wavenumbers  and  the  line  positions  by  less  than  a  wavenumber, 
are  insignificant  in  the  8  cm  ^  spectrally  resolved,  spectra. 
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The  computer  program  which  was  written  is  thus  severely  limited 
in  application.  It  can  only  reproduce  the  first  three 
vibrational  hot  bands  of  the  band  of  H2O  and  their  Q  branches. 
It  may  only  calculate  with  any  degree  of  certainty  low  resolution 
spectra  (8  cm~^) .  The  program  includes  m  provisions  for  higher 
perturbations  to  the  line  positions  or  line  intensities. 

The  program  used  the'  Scientific  Subroutine  "EIGEN"  for 
diagonalizing  the  energy  matrices.  The  program  calculates  the 
eigen  values  for  J  maximum  of  30. 

5.1.2  Slit  Function 

To  reproduce  theoretically  an  experimentally  obtained 
band  shape,  it  is  necessary  to  use  a  theoretical  instrument 
function  that  closely  resembles  the  actual  instrument  function. 

In  practice,  workers  have  successfully  used  both  Gausian  and 
triangular  slit  functions.  In  order  to  determine  which  function 
best  fits  our  instrument  function,  the  instrument  slit  function 
was  determined  experimentally.  It  was  observed.  Figure  24, 
that  a  triangular  function  more  closely  resembled  the  low 
resolution  instrument  function  than  did  the  Gausian  function. 

This  instrument  function  was  determined  by  scanning  a  laser  line 
of  width  =0.1  cm~^  with  the  monochromator  set  at  8  cm  There¬ 
fore,  the  triangular  slit  function  was  used  in  these  present 
calculations . 

5.1.3  Computer 

The  computer  program  was  written  in  Fortran  IV  language 
for  use  on  a  Data  General  -  NOVA  computer.  Since  the  core  size 
of  this  computer  is  small,  it  was  necessary  to  break  up  the  large 
program  into  four  smaller  ones.  The  NOVA  computer  with  disk 
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FIGURE 

CALCULATED  ‘RAMAN  BAND  PROFILES  OF  H^O 
IN  THE  REGION  OF  H^O  VAPOR  (8  CM~^  RESOLUTION) 
AT  VARIOUS  TEMPERATURES  RELATIVE  TO  298  K. 
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A.  CALCULATED  PROFILE 

B.  V(I00)-V(000)  COMPONENT 

•500.  TEMPERATURE  (K)  C.  V  (NO  )-V(OIO)  COMPONENT 

D.  V(I20)-V (020)  COMPONENT 


FIGURE  26 

THE  VIBRATIONAL  HOT  BAND  COMPONENTS  OF  H^O 
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5.4  THEORETICAL  LINE  ANALYSIS 

The  calculated  Raman  band  contours  of  water  vapor  at  1500  K 
are  shown  in  Figures  27  to  29  for  the  (100)-*- (000)  transition 
with  slit  functions  of  8,  4  and  2  cm  ^  respectively.  The  low 
resolution  contour  lacks  fine  structure,  however  at  increasingly 
higher  resolution  a  definite  fine  structure  which  is  only 
slightly  regular  appears. 

The  pattern  which  is  quite  clear  in  the  2  cm  ^  slit  function 
contour  is  unexpected  and  not  predicted  as  a  highly  asymmetric 
top  such  as  H2O  should  not  have  a  regular  pattern  in  its  Q  branch. 

To  determine  the  reason  for  the  peculiar  fine  structure  in 
the  2  cm”^  contour  of  Figure  29,  it  is  necessary  to  analyze  which 
lines  and  therefore  which  J  and  components  are  contributing 

to  the  fine  structure.  In  Table  10  are  the  component  lines 
which  are  found  to  be  contributing  to  the  intensity  of  peaks 
located  at  3631  and  3627  cm”^  of  Figure  29.  Note  that  the 
pattern  is  irregular  and  that  the  grouping  of  the  lines  to  form 
the  fine  structure  in  the  calculated  contours  is  merely  coinciden¬ 
tal. 

The  overall  intensity  of  the  band  contour  arises  from  a 
Boltzmann  distribution  within  the  Q  branch.  In  general,  it  is 
seen  that  those  lines  with  low  J  and  K  contribute  most  to  the 
band  intensity.  The  fall-off  in  intensity  as  J  or  K  increases 
is  strictly  a  statistical  result.  Furthermore,  the  shift  of  the 
line  from  the  band  center  increases  with  J  and  K.  At  temperatures 
above  1500  K  band  contour  extends  several  hundred  wave  numbers 
below  the  band  center.  An. analysis  of  the  line  positions  and 
intensities  indicated  that  very  little  intensity  was  contributed 
to  the  overall  band  intensity  for  J  values  greater  than  20  and 
Jmax  “  25  was  more  than  sufficient  for  calculating  the  band  con¬ 
tour  at  2000  K. 
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RAMAN  SPECTRUM  OF  H^O 
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FIGURE  28 

CALCULATED  RAMAN  PROFILE  OF  H^O ' 
WITH  4  CM"^  resolution 
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FIGURE  29 

CALCULATED  RAMAN  BAND  PROFILE 
OF  H»0  WITH  2  CM"^  RESOLUTION 
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TABLE  10 

PARTIAL  LIST  OF  LINE  POSITIONS 


Lines  in  Region  Lines  in  Region 

(3629.5-»-3632.5  cm"^)  (3627 . 5-»-3629 . 5  cm"^) 


J 

1 

J 

K.i 

7 

(0,7) 

(7,0) 

8 

(1,8) 

(5,4) 

8 

(1/7) 

(5,3) 

9 

(2,8) 

(6,3) 

9 

(2,7) 

(7,3) 

10 

(3,8) 

(8,3) 

10 

(3,7) 

(8,2)  (9,2) 

11 

(4,8) 

(9,2) 

11 

(4,7) 

(10,2) 

12 

(5,7) 

(12,1) 

12 

(6,6) 

(12,0) 

13 

(6,7) 

(13,0) 

1.  The  lines  in  the  3629. 5*>-3632.5  cm  ^  region  correspond  to  a 
peak  in  the  2  cm  ^  resolution  spectra  -  Figure  29. 

■2.  These  lines  in  the  3627.5->-3629. 5  cm  ^  region  correspond  to 
a  valley  in  between  peaks  in  Figure  29. 
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5.5  COMPARISON  OF  CALCULATED  CONTOURS  TO  OBSERVED  CONTOURS 

The  comparison  between  the  observed  Raman  band  contour 
and  the  calculated  contours  is  demonstrated  in  Figure  30.  The 
observed  Raman  band  profile,  obtained  at  -  1100  K,  is  plotted 
with  the  calculated  contours  of  1000,  1100  and  1200  K.  There 
appears  to  be  excellent  agreement  between  the  calculated  tempera¬ 
ture  profile  at  1100  K  and  the  temperature  obtained  from  the  H2O 
Raman  spectrum.  These  results  clearly  show  that  if  the  signal 
to  noise  statistical  error  is  sufficiently  low  the  flaune  tempera¬ 
tures  obtained  from  an  analysis  of  the  Raman  band  of  H2O  are  com¬ 
parable  to  those  previously  obtained  from  analysis  of  the  N2  band. 

The  experimental  spectra  discussed  in  Section  III  have  an 
uneven  band  contour  which  appears  to  increase  with  temperature 
and  at  first  examination  appears  to  be  fine  structure.  However, 
the  spectra  obtained  at  1650  and  1750  K  do  not  have  this  shape 
in  their  contours.  It  is  suspected  that  this  structure  on 
the  contoxir  arises  from  Na  emission:  the  spectra  obtained  be¬ 
tween  1000  and  1500  K  were  carried  out  by  photon  counting  tech¬ 
niques,  whereas  the  higher  temperature  spectra  were  obtained  using 
phase  sensitivity  detection  which  eliminated  the  flame  emission 
problem. 

The  experimental  Raman  band  profiles  show  quite  clearly 
the  relationship  between  the  asymmetry  and  the  temperature  of 
the  water  band.  A  relationship,  therefore,  should  exist  between 
the  full  width  at  half  height  (FWHH)  of  the  Raman  band  profile 
and  the  rotational-vibrational  temperature  of  the  water  molecules. 
In  Figure  31  a  theoretical  FWHfl  versus  temperature  plot  is  shown. 
The  experimentally  observed  FWHH's  of  the  water  band  at  the 
temperatures  determined  from  the  N2  band  analysis  have  been  put 
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EXPERIMENTAL  SPECTRUM 

H2O  (nOO  K)  CALCULATED  SPECTRA 
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-*-1100  K 
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FIGURE  30 

COMPARISON  OF  EXPERIMENTAL 
AND  CALCULATED  RAMAN  BAND  PROFILES 
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FIGURE  31 

PLOT  OF  TEMPERATURE  VERSUS  THE  FULL  WIDTH  AT  HALF  HEIGHT 


OF  THE  EXPERIMENTAL  AND  THEORETICAL  RAMAN  BAND  CONTOURS  OF  H^O  VAPOR 
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in  the  figure  to  show  the  comparison.  The  error  bars  shown  on 
this  data  represent  the  one  sigma  error  uncertainty  in  deter¬ 
mining  the  FWHH  (±25  K) .  Note  that  within  experimental  error, 
the  use  of  the  experimentally  obtained  FWHH's  is  accurate  enough 
for  equilibrium  temperature  determination  of  the  water  molecules 
in  the  flame  gases. 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 

‘6.0  SUMMARY  OF  WORK  ACCOMPLISHED 

In  accordance  with  the  Work  Statement  for  this  contract, 
the  work  accomplished  and  submitted  in  this  technical  report 
includes,  but  is  not  limited  to,  the  following: 

a.  Experimental  Study 

(i)  Results  have  included  Raman  Stokes  cross  sections 
of  H2O  relative  to  N2  at  room  temperature  and  at 

a  minimum  of  five  temperatures  between  1000  -  1750  K. 
Section  III. 

(ii)  Verification  that  temperatures  were  calibrated  us¬ 
ing  Raman  scattering  of  nitrogen.  Section  III. 

(iii)  Estimates  of  inhomogeneity  of  the  source  tempera¬ 
ture  and  the  flame  specie  ntimber  density  values 
are  discussed  in  Section  III.  Due  to  the  large 
number  of  variables  present  in  the  experiments, 

a  measurement  of  these  errors  was  too  difficult 
to  obtain. 

(iv)  Verification  of  the  wavelength  calibration  of 
the  monochromator  and  that  the  wavelength  uncer¬ 
tainty  is  inconsequential  to  the  spectra  obtained. 
Section  III.  The  instrument  was  calibrated  with 
lamp  standards  to  an  accuracy  of  ±1  cm~^  and  since 
the  present  data  was  low  resolution  ('''  8  cm  the 
wavelength  uncertainty  can  be  disregarded. 

(v)  The  major  contributions  to  the  experimental  error. 
Section  III. 

(vi)  Documentation  that  studies  were  performed  at  a 
second  wavelength  to  verify  that  Raman  cross-section 
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scaling  varies  as  the  source  wavelength  to  the 
fourth  power.  Section  III. 

b.  Theoretical  Study 

(i)  A  computer  program  (in  Fortran  IV  code)  is  pro¬ 
vided  that  allows  for  only  the  theoretical  calcu¬ 
lation  of  the  water  vapor  band  Q  branch  profile 
at  various  temperatures.  Appendix. 

(ii)  Calculated  Raman  band  contours  of  water  vapor  in 
the  region  over  a  wide  range  of  temperatures. 
Section  V. 

(iii)  Experimental  verification  that  the  contribution 
of  the  anisotropic  Raman  scattering  by  the 
band  of  water  vapor  is  negligible  and  need  not 
be  considered.  Section  III.  This  is  true  only 
for  these  low  resolution  studies.  Polarization 
experiments  at  high  resolution  are  required  for 
further  verification. 

(iv)  The  results  of  the  theoretical  study  include  an 
interpretation  and  correlation  of  the  water  vapor 
spectrumr  an  analysis  of  the  observed  Raman  band 
energy  levels  in  the  region  of  water  vapor,  and 
a  comparison  of  the  observed  experimental  Raman 
band  contours  and  the  calculated  contours. 

6.1  USE  OF  H2O  RAMAN  BAND  PROFILES  FOR  TEMPERATURE  AND  SPECIE 

DENSITY  STUDIES 

The  results  of  this  study  show  quite  clearly  that  Raman 
scattering  studies  of  water  vapor  can  give  temperature  and  specie 
densities.  The  work  has  demonstrated  that,  if  signal-to-noise 
statistics  allow,  the  Raman  data  obtained  from  water  vapor  gives 
comparable  results  to  that  obtained  from  N2  Raman  data.  The 
accuracy  of  these  determinations  is  not  greater  than  +  50  K  for 
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temperature  and  not  greater  than  ±10%  in  number  density. 

6.2  RECOMMENDATIONS  FOR  EXTENDED  WORK 

6.2.1  Modification  of  the  Computer  Program 

Since  this  work  clearly  demonstrates  the  high  degree  of 
success  in  reproducing  theoretically  the  observed  Raman  band 
contours  of  low  resolution  H2O  spectra,  it  is  recommended  that 
future  work  be  carried  out  to  extend  this  study.  The  study 
should  include  the  following: 

(1)  the  generalization  of  the  computer  program  for 
calculation  of  Raman  spectra  for  any  molecule  of  any 
symmetry  in  the  vapor  phase; 

(2)  a  theoretical  study  to  determine  the  intensity 
expressions  for  the  anisotropic  Raman  scattering  of 
asymmetric  tops; 

(3)  modification  of  the  program  for  calculation  of 
theoretical  high  resolution  spectraT  This  will  include 
the  ability  to  include  the  effects  of  Doppler  broadening, 
pressure  broadening,  centrifugal  distortion,  Fermi 
resonance,  Coriolis  coupling  effects,  and  several  kinds 

of  slit  functions  in  the  theoretically  calculated  spectra. 

6.2.2  Experimental  Work  on  Other  Molecules 

It  is  recommended  that  Raman  scattering  experimental 
studies  be  carried  out  on  other  flame  species  of  interest 
simultaneously  with  the  theoretical  studies  proposed.  These 
studies  should  be  carried  out  to  obtain  cross  sections  of  interest 
and  to  verify  the  theoretical  calculations.  The  possibility  of 
using  Raman  scattering  spectroscopy  in  burnt  fuel  diagnostics 
has  been  clearly  demonstrated  by  a  nxamber  of  workers  and  by  this 
study.  Advantage  should  be  taken  of  this  technique  and  it  should 
be  pursued  to  its  fullest. 
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APPENDIX ■ 

COMPUTER  PROGRAM  FOR  CALCULATING 
LOW  RESOLUTION  RAMAN  BAND 
PROFILES  OF  H2O  IN  THE 
REGION 
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(J 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


2 


THIS  bfute  force  pfogfam  calculates  the 
TFaCE  SCATTEFING  RAMAN  BAND  PROFILE  SPECTRUM 
FOF  H20  AT  VAPIOUS  TEMPERATURES  UP  TO  2e00K 

C 

DIMENSION  AC  1  326  > «  A 1  C 6  1 Aec  21 «  C0C  21«  ACC800> 

THE  FOLLOWING  SERIES  OF  DO  LOOPS  COMPUTES  THE  ELEMENTS 
OF  THE  ENERGY  MATRICES  FOR  THE  VIBRATIONAL  STATES  OF  H20 
C000>«  (  102)4  (0iei«  C  1  10)«  (e20WC  1201  FOP.  VALUES 
OF  J  UP  TO  30  ■ 


CALL  RESET 
CALL  FOPEN 
CALL  FOPEN 
CALL  FOPEN 
CALL  FOPEN 
CALL  FOPEN 
CALL  FOPEN 


C04 ' A000*/A1 
C  l4  '  A100*4  41 
(24  'A010S41 
C  34  'Al  10*44) 
(44  *  A020*«  4) 
(54  •A120*4  4) 


1  Ils0 

DO  10  M»  14  6 
A0(  11  =  27.795 
B0( 11= 14.508 
C0( 11=9.289 
A0C21  =  27. 300 
50(21=  14.284 
C0( 2)=9. 144 
A0(  31  =  30.  454 
B0( 31= 14. 170 
C0(3)=9.  184 
A0(4)  =  29.959 
50(41= 14.486 
C0(4l=9.039 
A0(5)  =  33.  1  1  3 
50(51= 14.9  12 
C0(5l  =  9.879 
A0C6)  =  32.618 
50(61= 14.688 
C0(61=8.934 
L=25 


LLs2 

DO  1  J=  I4L 

ME=((  (2*Jl-t>  l)*=2-^C(2=Jl-^  111/2 
DO  2  JJ=  I4ME 


A(JJ1=0.0 

CONTINUE 
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ME*  <  <  C  2*J  » +  I  >»*  2+C  (  2»J>  ♦  1 ) » /  2 

II*  I 

K*-J 

MEE*2»J+2 
DO  3  I*2/MEE 

ACIIJ  =  .5*CB0(M>  +  CISCMJ»*tJ*tJ+l>-tK»*2)  >  +  A0CM  J*CK**2» 
IFCK. EC. -  I >  A< 1 1 )*  A(  1 1 £5«C  B0CM>- C0CM> )*J«C J+ I J 
IFCK.EO.  1>ACII>  =  ACI1>  +  . 25* «  B0CM » - C0tM >  > *J*< J+  U 
IFCK.LT.0)GO  to  33 

'  ACll>=AClI)+0. 1 
33 

K*K+  1 

3  continue 
11*4 
K»-J 

DO  4  I*4>MEE 
N*J»C J+ 1 > 

ZT*  1. 0*CN-K*CK+ 1 > 1 
ZP.=  1. 0*CN-CK+  l)*CK+2)  ) 

AtIIi=.25»tE0«M>“C0tMi i*tZT»ZP)**. 5 

IF  CK.EC.-2J  AC1IJS0.0 

IF  CK.EC.-U  ACII>  =  e.0 

IFCK. EG. 0J  ACI I J*C2.0»*.5J*AC1I  > 

1 1=11  +  1 
K*K+  1 

4  CONTINUE 
C 

c  the  FOLLOVING  eefies  of  statements  calls 

c  EUBPOUTINE  EIGEN  PNC  DIAGONALIZES  THE 

C  ENEPGY  MATRICES  FOP  THE  ENEFGY  LEVELS  OF 

C  E/\,CH  VIBRATIONAL  STATE  AND  DETEFMINES  THE 

C  TPANSl.TlONAL  ENERGIES  FOR  THE  0  BRAUCH  OF  THE 

C  FUNDAMENTAL  AND  THE  FIRST  TWO  HOT  BANDS 

C 
C 

N*J 

IFCJ.EG.il  GO  TO  16 
CALL  EIGENCA^R^N#  1> 

I  I*  I 

DO  5  I*  1> J 
AIC1>*ACII> 

II*II+I+1 

5  CONTINUE 

16  AlC  1  i*AC  1  1 

IFH*CJ*J+J)/2 
NN*J+IFH+1 
NU*J+  1 
MM*  1 

DO  11  I=1/NU 
DO  12  II*  1j  I 
ACCMM)*ACNN1 
MM*MM+ 1 
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NN»NN+ 1 

12  CONTINUE 
NN*NN+J 

11  CONTINUE 

NeJ+  I 

DO  14  1=  1<NU 
14  A(I)bACCI) 

CALL  EIGEN  CA^P^N^I* 

MMe  1 
JJ»2 

N»  C  2* J ) ♦  1 
DO  13  I  =  NU-N 
AlCI  )  =  ACMM> 

MM=MM+JJ 

13  JJ=JJ+1 

CALL  FSEEK  <1II^LL> 

DO  6  1-  1»N 
A2»A1CI  ) 

VFITE  BINAPV  CIIU  A2 
LL=LL+l 

6  CONTINUE 

1  CONTINUE 

CALL  FCLOS  till* 

IIIsI  II-»  1 
TYPE  'LEVEL ’#M 

10  continue 

CALL  FCHAN  (’TEST2.SV'> 
END 
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SUBROUTINE  EIGENC  A/ R^NjMV) 

CIMENSION  AC  1  PC  1 ) 

5  FANGE®1.0E-6 

IFCMV-  1)  10>  25>  10 

10  IG=-N 

DO  20  UN 
I  Qs  I  Q+N 
DO  20  1»  UN 
IJBl  C-Kl 
P.CIJ)B0.0 
IFCI'J)  20*  15/20 
15  P.ClJ)sU0 

20  CONTINUE 

25  /^OPM::0*0 

DO  35  I»  UN 
DO  35  J=I/N 
IFCI-J}  30/  35/  30 
30  I A»I+tJ«J-JJ /2 

ANOPM»ANOPM  +  ACI  A>*ACIA» 

35  CONTINUE 

1FCAN0PM>  165/  165/40 
40  AN0F-M=  U414«SGRTCAN0RM) 

AN  FMXb  ANOPM*  RANGE/ FLO  ATCN  » 

INDS0 

THR=AN0RM 

45  THF«THR/FL0ATCNJ 

50  L»  I 

55  M»L+ 1 

60  MG=CM*M-M>/2 

LQ=CL*L-L)/2 
LM=L+MQ 

62  IFC  ABSCACLM  J  J-THRJ  130/65/65 

65  INI>1 

LL-L+LQ 

MM=M+MG 

Xs0.5«CACLL)-ACMM) ) 

68  Ya-ACLM>/  SCRTCACLMX'ACLMifX-Xi 

IFCXl  70/75/75 
70  Y=-Y 

75  £INX=Y/  £QRTC2.0*C 1.0+C  SGRTC 1 . 0-Y*Y > > > > 

?INX2=  £1NX*S1NX 
78  COSXa  SORTC  1.0-S1NX2> 

COSX2=COSX*COSX 
SINCS=£1NX*C0£X 
ILO=N*CL-  1  J 
IMQ=N"«CM-  1> 
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DO  125 

lG=Cl*I-I>/2 
IFCI-L>  80<  1  15^82 
80  85*  1  15*90 

85  IM=I+MQ 

GO  TO  95 
90  IM=M+IG 

95  1F<1-L}  100*  105* 105 

100  1L»1+LQ 

GO  TO  110 
105  1L=L+IQ 

110  X=AtIL)»COSX-ACIM)*SINX 

ACIM)  =  A«  ILJ*SIKX  +  AC1M)*C0SX 
A(IL)sX 

1  15  IFCMV- 1 J  120*  125*  1 20 

120  ILF=ILe+i 

lMF=lMn+I 

X=P.(  ILFi*CO£X-FCIMR)*SINX 
PC  IMF)  =  PCILR)*£INX  +  ECIMF.J*COSX 
FC1LP>=X 
125  CONTINUE 

X=2.0*ACLM)»S1NCS 
Y=ACLLJ*C0£X2+ACMM>*£INX2-.X 
X=AtLLJ*SINX2+ACMMJ»CO  £X£+X 

ACLM)=  CACLLJ-ACMMJ  1 »  SINC£+ ACLM  C  CO  £X2- SINX 2  J 

ACLLJ=Y 

ACMM>»X 

130  IFCM-Ni  135*  U0*  135 

135  M=M-M 

GO  TO  60 

U0  IFCL-<N-1))  145*150*145 

145  L=L+I 

GO  TO  55 

150  IFCIND-1>  160*  155*  160 

155  1ND=0 

GO  TO  50 

160  IFCTHP-ANRMX)  165*165*45 

165  IC=-N 

CO  185  I=1*N 
IC=IG+N 

LL=I+t 1*1-1 >/2 
JCi=N»tI-2J 
DO  185  J=1*N 
JC=JC+N 

MM=J+C J*J-J>/2 
IF<  AtLLJ-A<MM»  >  170*  185*  185 

170  X=ACLLJ 

A<LL»  =  A<MM» 

A«MM)=X 

IFCMV-  1  >  175*  185*  175 

175  CO  180  K= 1*N 
lLr=IC+K 
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IMF:^JC+K 
X«F<1LF.) 
FCILF>  =  F.C1MPi 
PCIMRl=X 

continue 

PETUPN 
END 
F. 
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DIMENSION  F1NTC400> 
COMMON  FINT,  TEM 
TYPE  ‘TESTg* 

CALL  RESET 

CALL  FOPENC04  'Aee0*«  4) 
CALL  FOPENC  1^  'A100'«  4) 
CALL  FOPENC  'AO  10*«  41 
CALL  FOPENC  3#  '  A1  10S41 
CALL  FOPEN  C4#*A0S0*«4> 
CALL  FOPEN  (5#'A120'«41 
L=25 

LN=L+ 1+CL>»CL+ 11 
EL1=3657.054 
ELI  1»EL1 

EL2S5225*  0-  1590.0 
EL21  =  EL2 

EL 3=676  1. 0-3151.0 
EL31  =  EL3 
DO  1  1  IJ=2>LN 
CALL  FSEEK  C 1«  XJl 
READ  BINARY  <11  A2 
A22=A2 

CALL  FSEEK  C0^IJ1 
READ  binary  <01  A2 
EL1=A22-A2+EL  1  1 
CALL  FSEEK  <  1«  IJl 
WRITE  BINARY  <13  ELI 
CALL  FSEEK  <3.  IJl 
READ  BINARY  <3l  A2 
A23=A2 

CALL  FSEEK  <2,  IJl 
READ  BINARY  <2l  A2 
EL2=A23-A2+EL2  1 
CALL  FSEEK  <3#  IJl 
WRITE  BINARY  <31  EL 2 
CALL  FSEEK  <5«IJ1 
read  binary  <51  A2 
A24= A2 

CALL  FSEEK  <4^ IJl 
READ  BINARY  <41  A2 
EL3=A24-A2+EL31 
CALL  FSEEK  <5#1J1 
WRITE  BINARY  <51  EL3 
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continue 

CALL  FSEEK  C  1j  11 
WRITE  BINARY  <11  EL  1  I 
CALL  FSEEK  C  3j  1 1 
WRITE  BINARY  <31  EL21 
CALL  FSEEK  C 5# 1 1 
WRITE  BINARY  <51  EL 31 

the  REMAINER  of  THE  PROGRAM  CALCULATES  THE 
LINE  INTENSITIES  AND  CONVOLUTES  THE  LINE 
SPECTRUM  INTO  AN  INSTFIMENT  SPECTRUM 
WITH  A  TRIANGLI-AF  SLIT  FUNCTION  USING  A 
BANDWITH  OF  8  WAVE-NUMBERS  FOR  VARIOUS 
TEMPERATURES 

ICJIIS  PROPORTIONAL  TO  SCJ+HG  EXPC- C  EROT+EVI  Bl /KT) 

CONST=C6.62559E-27*‘3.0EI01/C  1. 3805AE-  16> 

TEM- 1500. 0 
ELI  0.0 
EL21sI590.e 
EL3I»315I.0 

EXPl  1  =  EXPC-(EL1  ll*CON£T/TEMl 
EXP2 1  =  EXPC -  t  EL 2  1 1  •  CON  ST/TEM  > 

EXP3  1  =  EXP<  -  «  EL  3  1 1  *  CON  ST/TEM  1 
DO  12  IJs2«LN 
CALL  FSEEK  C0^IJ) 

READ  BINARY  C0>  A2 
A22=A2 

CALL  FSEEK  C2«IJ] 

READ  BINARY  C2)  A2 
A23=A2 

CALL  FSEEK  CA^IJl 
READ  BINARY  CAl  A2 
A2A=a2 

EXPl=EXPt-CA22+EL  1  1  >»  CON  ST/TEM  > 

CALL  FSEEK  C0«IJ1 

WRITE  BINARY  t0>  EXPl 

EXP2»  EXPC-«A23+EL21)*  CON  ST/  TEM  ) 

CALL  FSEEK  C2/1J1 
WRITE  BINARY  <2)  EXP2 
EXP3=EXPC-C  A2A+EL3I)*C0NST/TEM> 

CALL  FSEEK  C4«IJ) 

WRITE  BINARY  £4l)  EXP3 

CONTINUE 

CALL  FSEEK  C0«l} 

WRITE  BINARY  C0>  EXPl  I 
CALL  FSEEK  C 2« I 1 
WRITE  BINARY  C2)  EXP21 
CALL  FSEEK  C  4^  1 1 
WRITE  BINARY  C4)  EXP3i 
N«0 
K=  1 
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LI=L+ 1 

CALL  FSEEK  <0*  1> 

CALL  FSEEK  C  2,  1 1 
CALL  FSEEK  C4«l} 

DO  1  3  tJ»  l^LI 
JJ=<2»J)- I 
DO  22 

read  binary  C0J  EXPl 
EXPJ  I  =  EXP1#CC2*NJ+  11 
CALL  FSEEK  C0,K> 

VP.ITE  BINARY  C0>  EXPJ  1 
READ  BINARY  C2l  EXP2 
EXPJ2=EXP2*CC2»Nl  +  1  > 
CALL  FSEEK  C2,K) 

WRITE  BINARY  <21  EXPJ 2 
read  binary  C41  EXP3 
EXPJ3=EXF3*(  (  2*N>  +  1  > 
CALL  FSEEK  C4«K> 

WRITE  BINARY  C4)  EXPJ  3 
K=K+1 

22  CONTINUE 
N=N+J 

13  CONTINUE 
LI=L+1 

CALL  FSEEK  C  0,  1 1 
CALL  FSEEK  C  2#  1 1 
CALL  FSEEK  (4,1) 

K=  1 
G=3.0 
G  1=  1.0 

DO  14  J=  1,LI 
JJ=C2»J)-  1 
DO  23  1=1# JJ  ■ 

READ  BINARY  (0)  EXPJ  I 

tensi=expj  i#g 

CALL  FSEEK  (0#K} 

WRITE  BINARY  C0)  TEN  S  1 
READ  BINARY  (21  EXPJ 2 
TENS2=EXPJ2*G1 
CALL  FSEEK  (2#K) 

WRITE  BINARY  (2)  TENS2 
READ  BINAP?Y  (4)  EXPJ  3 
ten  S 3=  EXPJ  3* G 
CALL  FSEEK  (4#K1 
WRITE  binary  (41  TENE3 
J  1=J-  1 

IFd.NE.Jll  GO  TO  230 
GG=G 
G=G1 
G  I  =  GG 
230  K=K+1 

23  CONTINUE 

14  CONTINUE 

CALL  FCHAN  ( 'TESTS.  SV'l 
END 


F 
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CONVOLUTION  STEP  FOP  3500-3702  VAVE-NIMBEP.S  PEGION 

DIMENSION  FINTCASei  #  AEL  K  7001  /  AEL2C  700  AEL  3<  700  ) 
DlMeJSION  ATENSIC700J  /  ATEN'S2C  7001 ,  ATENS3C  700 J 
COMMON  FINT/TEM 
TYPE  'TEST3’ 

CALL  FOPEN  C0,  'A000'/A) 

CALL  FOPEN  CU*A100’*4i> 

CALL  FOPEN  (2^’A010'#A) 

CALL  FOPEN  C3^'A110’#4) 

CALL  FOPEN  (4,  •A020SA> 

CALL  FOPEN  <5j’A120S4) 

L=25 

LN=L+  1  +  CL>«CL+ 13 
CALL  FSEEK  C0,  13 
CALL  FSEEK  C  1,  13 
CALL  FSEEK  C  2/  1  3 
CALL  FSEEK  C  3j  1  3 
CALL  FSEEK  t 4,  1  3 
CALL  FSEEK  C5^  13 
DO  1  1=I^LN 
READ  BINAFY  C03  TEN  S 1 
READ  BINARY  C  13  EL  1 
READ  BINARY  C  23  TEN  S2 
READ  BINARY  £33  EL 2 
READ  binary  £43  TEN  S3 
READ  BINAFY  £53  EL 3 
AELU13SEL1 
AEL2£I  3  =  EL2 
AEL3£I  3SEL3 
ATENSUI  3  =  TEN  SI 
ATENS2£I  3  =  TEN£2 
ATENS3£  I  3  =  TENS3 
1  CONTINUE 

CONST*  1.43978  22 

AP.EA=0.0 

DO  40  I*  ULN 

40  AREA=APEA+ATENS1CI  3+aTENS2£I3+aTENS3£13 
DO  41  I=1>LN 

ATEN  S 1  £  I  3  =  ATEN  S 1  £  I  3 /AREA 
ATEN£2£I  3  =  ATENS2£I  3 /AREA 

41  ATENS3£I  3*aTENS3£I  3/aREA 
lNC=45e 

DO  23  1=  1#  INC 
23  FINT£I3  =  0.0 

MIN=350e 
MAX*  3700 
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SVsS. 

FW=  2* SU 
DO  25  K=1*LN 
KEL  1  =  AEL  1CK> 

KEL2=  AEL2CK) 

KEL3=AEL3CK1 
I=«KEL1-MIN)»2 
I I=CKEL2-MIN)*2 
III  =  CKEL3-MINJ*2 

I FC I . GT. 4 I6.0R. I •LT.- I6> GO  TO  28 
LM=-  16 

DO  26  ILM=LM>16 

WT 1  =  ATENS  IC  K>  *  C  FV?-  I  ABSt  ILM  J  >  /FV 
IFCWTl.LE.e.0>GO  TO  26 
M=I+ILM 

lFtM.LE.e)GO  TO  26 
FINTCM)  =  FINTCM)+WT1 

26  CONTINUE 

28  IFC  1 1  .GT.  4  16.0R.  1 1  .LT.- 16>60  TO  29 
LM=-I6 

DO  27  ILM=LM^ 16 

VT2=ATENS2CKJ»CFU-IABS<1LM>  >/FU 
IFCWT2.LE.0. 0>GO  TO  27 
M=I1+ILM 

IFCM.LE.0JGO  TO  27 
F1NTCM>=FINTCM)+WT2 

27  CONTINUE 

29  IFCIII.GT.A16. OF. III. LT.eiGO  TO  25 
LM«-  16 

DO  30  ILM=LM>  16 

WT3=ATENS3CK1*CFW-1 AESCILMJ  >/FW 
IF(VT3.LE.B.  0>GO  TO  30 
M=1II+1LM 

IFCM.LE. 0. 0) GO  TO  30 
FINT<M)  =  FINTCM>+WT3 

30  CONTINUE 

25  CONTINUE 

ANORM=. 612553 
DO  A01  Is  I,  400 

401  FINTCI  IsFINTCl  >/aNOPM 

CALL  FCHAN  C*TEST4.£VM 
E30D 
R 
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C  riOFILE  FOr  A  SPECIFIC  tempefatufe 

c 

DIMENSION  LEL  1(8j*LEL2<  1 0  >  ,  IHOL  1 9  >  #  LHOL  C  12> 
DIMENSION  FINT  C45e> 

COMMON/LE/LBL  I*  LEL I HOL*  LHOL 

COMMON  FINT,  TEM 

data  LELl/'FFECUENCY  CM- I  V 

data  LEL2/ •  RELATl  VE  INTENSITY*/ 

data  IHOL/ 'TEMPEFATUFE  CKi  */ 

data  LHOL/'PAMAN  SFECTF.IW  OF  H20  */ 

TYPE  'TESTA* 

YY=  0. 0 
XM=  e . 0 
FLTXa  3500.0 
FLTY=e.0 

CALL  INITAL  (6>  200^  1  1.0/  17.0) 

CALL  PLOT  C0.5,  l.2,-3> 

CALL  PLOT  <0.0, 0.0,0) 

CALL  AXIS  t  0. 0,  0. 0,  LBL  1,  -  1  A,  1  0.  0,  0.  0,  3520. 
CALL  AXIS  C  0. 0,  0. 0,LBL2,  18,  5.  0,  90. 0,  0.  0,  0. 
CALL  Nl'MBEP  C  0. 5,  6.  0,  .  2  1,  TEM,  0.  0,  0> 

CALL  SYMBOL  C  A.  0,  6 . 0,  .  2  1,  IHOL,  0 . 0,  1  5) 

CALL  SYMBOL  C  0. 5,  7. 0,  •  2  l,LHOL,  0. 0,  2£) 

M=0.0 

CALL  PLOT  10.0,0.0,3) 

CALL  FENDN 
DO  73  1=  1,  A00 
Y=FINT< I )*5.0 
CALL  PLOT  CX,Y,  I) 

73  X=X+0.025 

CALL  FEMUP 
CALL  PSTP 
CALL  FESET 
100  STOP 

END 

r 
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